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Chapter  1 


Introduction 


The  rapid  development  of  the  state-of-the-art  in  monolithic  microwave  integrated  circuits 
has  intensified  the  need  to  develop  sophisticated  CAD  tools  for  use  in  device/circuit  analysis 
and  design.  In  particular,  to  increase  the  accuracy  and  effectiveness  of  CAD  tools  there 
is  a  need  to  develop  device  models  based  upon  fundamental  device  physics.  Active  device 
simulators  have  improved  rapidly  in  the  last  few  years  and  now  are  being  developed  to  the 
point  where  both  accuracy  and  simulation  efficiency  are  sufficient  for  practical  ffirst  pass’ 
design  applications,  that  is  device /circuit  designs  which  initially  prove  successful. 

In  order  to  obtain  the  maximum  benefit  from  a  device  simulator,  the  device  model  should 
be  capable  of  describing  the  performance  of  a  device  before  fabrication.  In  this  manner 
the  link  between  the  process  and  RF  simulators  would  be  closed,  and  much  time,  effort, 
and  expense  would  be  saved  since  device  optimization  studies  could  be  performed  before 
the  device  is  actually  fabricated.  This  consideration  indicates  a  physics  based  model,  and 
the  need  to  simulate  RF  operation  indicates  an  efficient  device  algorithm  based  upon  either 
smalytic  or  table  look-up  techniques.  Although  a  variety  of  physics  based  device  models  have 
been  reported  in  the  literature,  most  of  the  models  in  practical  use  are  based  upon  equivalent 
circuit  techniques  (e.g.,  the  MESFET  models  in  TOUCHSTONE  and  SUPER-COMPACT) 
and  require  that  the  device  be  fabricated  and  characterized  before  the  equivalent  circuit 
can  be  established.  This  procedure  does  not  allow  process  data  to  be  used  in  a  direct 
manner.  Also,  the  proper  definition  of  the  equivalent  circuit  requires  that  a  large  number  of 
devices  be  designed  and  fabricated  in  order  to  establish  satisfactory  and  sufficient  data  bases. 
The  accuracy  of  the  equivalent  circuit  techniques  is  not  well  established,  and  much  work 
remains  to  be  done  in  developing  suitable  device  cheuracterization  and  parameter  extraction 
techniques. 

A  physics  based  large-signal  MESFET  simulator  suitable  for  use  in  microwave  CAD  has 
been  developed  at  NCSU.  With  MIMIC  Program,  Phase  3,  support  the  simulator  has  been 
significantly  advanced  and  enhanced  with  an  RF  performance  and  process  yield  algorithm 
that  allows  the  simulator  to  be  used  for  both  device  design  and  process  yield  optimization. 
The  resulting  simulator  is  state-of-the-art  and  allows,  for  the  first  time,  the  physical  design 
of  MESFETs  to  be  optimized  so  that  a  maximized  percentage  of  devices  wiU  pass  a  spec¬ 
ified  yield  criterion  based  upon  large-signal  performance  such  as  RF  output  power,  gain, 
power-added  efficiency,  etc.  The  yield  simulator  defines  RF  performance  in  terms  of  the 
physical  design  parameters  actually  accessible  in  the  manufacturing  process.  The  simulator 
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has  proved  to  be  both  quantitatively  as  well  as  qualitatively  accurate  for  use  in  MIMIC  sim¬ 
ulations,  including  both  discrete  MESFET  devices  and  complete  integrated  power  amplifier 
circuits. 

The  development  of  the  simulator  has  been  facilitated  by  interaction  with  Phase  1, 
MIMIC  Program  participants.  In  particular,  experimental  data  has  been  supplied  by  ITT, 
Raytheon,  Texas  Instruments,  General  Electric,  and  Hughes.  Access  to  this  data  has  proved 
extremely  valuable  in  guiding  the  development  of  the  simulator  so  that  it  accurately  pre¬ 
dicts  device  performance.  The  simulator  has  been  verified  for  accuracy  against  a  variety  of 
industrial  devices,  with  varying  channel  doping  and  structural  designs,  operating  under  both 
clsiss  A  and  class  B  operating  conditions.  Examples  of  the  verification  tests  are  presented  in 
Chapter  2  of  this  report. 

Copies  of  early  versions  of  the  simulator  and  User’s  Manual  have  been  delivered  to  MIMIC 
Program  participants,  including  Raytheon,  Texas  Instruments,  COMPACT  Engineering, 
ITT,  Martin-Marietta,  General  Electric,  Hughes,  TRW,  Wright-Patterson  AFB,  and  Fort 
Monmouth.  The  simulator  and  manual  have  been  provided  at  no  cost.  In  addition,  user 
training  has  been  provided  to  ITT  and  Raytheon. 

The  device  model  and  verification  are  described  in  Chapter  2  of  this  report.  The  algo¬ 
rithms  for  computing  large-signal  performance  measures  and  using  them  for  sensitivity  anal¬ 
ysis,  nominal  performance  optimization,  and  3rield  optimizization  are  presented  in  Chapter  3. 
In  Chapter  4,  the  mathematical  optimizer  is  detailed.  Yield  optimization  experiments  are 
given  in  Chapter  5,  and  Chapter  6  closes  with  conclusions. 
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Chapter  2 

The  Physical  MESFET  Model 


The  simulation  philosophy  guiding  the  NCSU  simulator  is  indicated  in  the  flow  chart  for 
an  ’ideal’  and  integrated  CAD  procedure,  as  shown  in  Figure  2.1.  As  indicated  in  this 
chart  a  variety  of  models  are  required.  The  most  significant  factor  is  the  desirability  of 
interconnecting  process,  device,  and  circuit  simulators  so  that  a  MIMIC  can  be  simulated 
from  process  data  (such  as  device  geometry  and  channel  doping  details)  all  the  way  through 
to  RF  circuit  performance.  The  availability  of  such  an  integrated  simulator  would  allow,  for 
example,  a  device  design  to  be  optimized  based  upon  a  desired  RF  performance  specification 
(such  as  RF  output  power,  gain,  power-added  efficiency,  etc.). 

Various  building  blocks  for  the  ’ideal’  procedure  shown  in  Figure  2.1  are  under  develop¬ 
ment.  A  suitable  process  simulator  (e.g.,  SUPREM  3.5)  has  been  reported  [1].  A  version 
of  SUPREM  has  been  obtained  from  Stanford  University  and  the  code  has  been  linked  to 
the  NCSU  simulator  (TEFLON).  For  ion-implanted  devices  the  link  between  SUPREM  and 
TEFLON  permits  channel  doping  details  to  be  specified  directly  in  terms  of  process  param¬ 
eters  such  as  implant  species,  energy,  dose,  cap  thickness,  anneal  time,  etc. 

Suitable  microwave/mm- wave  linear  circuit  simulators  (i.e.,  the  integrated  circuit  models 
indicated  in  Figure  2.1)  have  been  intensively  developed  and  are  currently  in  a  relatively 
advanced  state  (e.g.,  the  commercially  available  simulators  TOUCHSTONE  and  SUPER¬ 
COMPACT).  These  simulators  also  include  versions  of  the  passive  element  models  indicated 
in  Figure  2.1.  The  linear  circuit  models  generally  work  fiurly  well  and  are  extensively  used 
by  industry.  In  TEFLON,  commercial  simulators  are  used  to  determine  characteristics  for 
necessary  passive  elements. 

2.1  Physical  Device  Models 

Accurate  device  models  suitable  for  use  in  microwave  CAD  simulators  can  be  formulated 
from  knowledge  of  the  device’s  physical  operation.  These  models  are  derived  from  solutions 
to  a  set  of  the  basic  semiconductor  device  equations  [2,  3],  subject  to  appropriate  boundary 
conditions  determined  by  the  MESFET  geometry  and  doping  profile.  The  main  difficulty  in 
applying  physical  device  models  to  microwave  CAD  simulators  is  the  trade-off  between  accu¬ 
racy  and  the  large  execution  times  generally  required.  Most  of  the  physical  models  solve  the 
semiconductor  device  equations  using  some  form  of  numerical  technique  such  as  finite  dif¬ 
ferences  [4,  5]  or  finite  elements  [6].  The  models  can  be  formulated  with  varying  complexity. 
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Figure  2.1:  Components  of  an  ‘ideal’  microwave/mm-wave  integrated  circuit  computer  aided 
design  system. 
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depending  upon  the  operational  detail  included.  For  example,  the  drift -diffusion  approxima¬ 
tion  [5]  where  the  charge  carriers  are  assumed  to  be  in  thermal  equilibrium  with  the  crystaJ 
lattice  is  often  assumed  for  devices  with  gate  lengths  on  the  order  of  a  micron  and  greater 
and  operating  through  X-band.  This  type  of  model  has  proved  very  valuable  in  illuminating 
operational  physical  phenomena.  More  advanced  formulations  include  energy-momentum 
relaxation  phenomena  to  account  for  nonequilibrium  effects  that  become  significant  for  sub¬ 
micron  gate  dimensions  and  very  high  frequency  operation  [7].  In  these  formulations  the  free 
electrons  are  not  in  thermal  equilibrium  with  the  crystal  lattice  over  at  least  a  portion  of  the 
conducting  channel.  Nonequilibrium  effects  tend  to  be  most  important  at  low  bias  where  the 
majority  of  the  carriers  are  in  the  low  effective  mass  conduction  band  central  valley.  As  bias 
or  RF  terminal  voltage  is  increased  the  carriers  transfer  to  higher  conduction  band  valleys 
where  higher  effective  mass  tends  to  dampen  the  hot  electron  behavior.  As  a  result,  hot 
electron  phenomena  is  reduced  for  power  devices  operated  in  saturation  where  the  terminal 
voltages  have  large  magnitudes. 

Quantum  effects  can  be  included  for  devices  such  as  HEMTs  by  including  Schrodinger’s 
equation  in  the  set  of  semiconductor  equations  [8].  Solutions  of  this  type  tend  to  be  complex 
and  require  extensive  computer  time  to  obtain  solutions.  The  models  have  not  yet  been 
developed  to  the  point  where  they  can  be  efficiently  used  to  simulate  the  RF  performance 
of  a  device.  In  general,  nonequilibrium  phenomena  can  be  simulated  using  hydrodynamic 
[9],  numerical  [7],  or  Monte  Carlo  [10]  solution  techniques.  These  models  are  very  useful  for 
investigating  in  precise  detail  the  physical  operation  of  the  device.  However,  simulation  time 
increases  rapidly  with  model  complexity  and  practical  operation  of  the  models  is  usually 
limited  to  DC  solutions.  Large-signal  operation  such  as  transient  performance  is  generally 
investigated  by  performing  a  series  of  DC  simulations  for  varying  electrode  potential.  There 
have  been  a  few  attempts  at  using  the  numerical  formulations  to  investigate  the  RF  operation 
of  a  device  [11, 12].  Conceptually,  there  is  no  reason  that  these  models  could  not  be  applied  to 
RF  simulations,  given  sufficient  computer  resources.  Such  simulations  are,  however,  beyond 
the  current  state-of-the-art. 

Hybrid  approaches  in  which  a  simplified  version  of  the  numerical  formulation  is  cou¬ 
pled  with  an  equivalent  circuit  approach  have  been  presented  [13,  14].  These  models  offer 
improved  accuracy  compared  to  the  simple  equivalent  circuit  models  and  provide  a  link 
between  the  numerical  simulations  and  RF  circuit  simulators.  The  basic  model,  however, 
remains  based  upon  numerical  techniques  which,  in  practice,  limits  the  designer’s  flexibility 
to  perform  device  and  circuit  optimization  studies  where  large  numbers  of  simulation  runs 
are  required.  For  example,  every  time  the  RF  circuit  is  tuned  a  complete  re-calculation  is  re¬ 
quired.  Since  many  tuning  conditions  are  required  to  determine  optimum  circuit  conditions 
simulation  time  quickly  becomes  unmanageable. 

Solution  of  the  device  equations  completely  by  analytic  methods  provides  a  good  com¬ 
promise  between  model  accuracy  and  execution  time  [15,  16].  In  this  approach  equivalent 
circuits  are  never  used  since  the  basic  model,  due  to  efficient  analytic  formulation,  can  be 
solved  directly.  When  this  approach  is  coupled  with  table  look-up  techniques  execution  time 
can  be  on  the  same  order  of  magnitude  or  less  than  that  required  for  the  equivalent  circuit 
models.  Since  device  nonlinearities  are  inherent  in  the  basic  semiconductor  equations  it  is 
not  necessary  to  make  a  priori  assumptions  as  to  the  form  or  identity  of  model  nonlinear 
functions.  The  basis  for  the  model  is  the  device  geometry,  channel  doping  details,  and  bias 
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and  RF  operating  conditions.  The  model  is  well  suited  for  both  device  and  RF  circuit 
optimization  studies.  Application  of  the  technique  to  RF  performance  and  process  yield 
optimization  is  attractive  since  large  numbers  of  simulations  can  be  efficiently  performed. 


2.2  The  NCSU  Microwave  MESFET  Simulator 

(TEFLON) 

The  NCSU  physical  GaAs  MESFET  simulator  (TEFLON)  is  based  upon  efficient,  ana¬ 
lytic  solution  of  the  basic  semiconductor  device  equations.  The  device  model  is  quasi-static 
and  is  based  upon  an  analytic  formulation  of  the  charge  dipole  domain  [15, 16].  This  formula¬ 
tion  has,  in  turn,  resulted  in  the  derivation  of  a  self-consistent  large-signal  analytic  MESFET 
model  that  can  easily  be  implemented  in  microwave  CAD  simulators.  The  time  domain  ana¬ 
lytic  device  model  is  interfaced  with  an  RF  circuit  by  means  of  the  harmonic  balance  method 
to  produce  an  efficient  and  accurate  nonlinear  device/circuit  microwave  simulator. 

The  model  permits  the  RF  performance  of  a  device  or  integrated  circuit  to  be  deter¬ 
mined  as  a  function  of  process  and  device  design  information  and/or  bias  and  RF  operating 
conditions.  In  this  manner  a  complete  integrated  device/circuit  simulator,  as  indicated  in 
Figure  2.1,  is  possible.  Such  an  integrated  simulator  allows  both  the  active  device  and  the 
passive  elements  to  be  optimized,  based  upon  the  parameters  accessible  in  the  fabrication 
process.  That  is,  factors  such  as  device  geometry,  ion-implant  species,  dose,  and  energy 
that  result  in  optimized  RF  output  power,  power-added  efficiency,  etc.  can  be  determined. 
This  approach  eliminates  the  need  to  experimentally  characterize  the  devices  in  order  to 
parameter  extract  equivalent  circuit  element  values.  The  experimental  measurements  that 
are  of  interest  consist  of  RF  performance  data  such  as  Pout  vs.  Pin,  PAE,  gain,  spectrum, 
and  impedance.  This  data  can  be  directly  compared  to  simulated  data  to  verify  the  model. 
Once  verified,  the  model  can  then  be  used  in  the  design  process. 

A  block  diagram  for  the  NCSU  large-signal,  analytic  model  indicating  the  information 
flow  is  shown  in  Figure  2.2.  The  inputs  can  be  supplied  from  experimental  measurements  or 
calculated  data.  The  outputs  of  the  simulator  consist  of  DC  I-V  characteristics,  small-signal 
parameters,  and  large-signal  data.  A  major  advantage  of  this  type  of  simulator  is  that  the 
DC  and  RF  performance  of  the  device  can  be  determined  as  a  function  of  parameters  that 
are  actually  accessible  in  the  fabrication  process. 


2.3  Verification  Experiments 

In  practice,  the  simulator  is  currently  limited  to  single  device  circuits  operating  in  the 
common-source  configuration,  as  shown  in  Figure  2.3.  Any  arrangement  of  source  and  load 
impedances,  and  parallel  or  series  feedback  networks  can  be  accommodated.  The  common- 
source  limitation  is  not  fundamental  and  was  implemented  for  convenience.  The  simulator 
could  be  generalized  with  slight  modifications  to  the  code.  However,  this  has  not  been  done 
due  to  time  restrictions. 

The  simulator  requires  that  channel  doping  details  be  known.  The  data  can  be  calculated 
from  process  models  (e.g.,  SUPREM  3.5  [1])  or  determined  from  measurements.  In  order 
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to  facilitate  determination  of  the  doping  profile,  a  parameter  extraction  technique  has  been 
developed.  The  parameter  extraction  technique  automatically  adjusts,  by  means  of  a  mathe¬ 
matical  optimization  procedure,  the  magnitude  of  the  physical  parameters  in  the  model  until 
good  agreement  between  a  variety  of  performance  characteristics  are  obtained.  Typically, 
these  consist  of  the  drain  and  gate  currents  in  addition  to  RF  output  power,  PAE,  and  gain. 
The  parameter  extraction  routine  has  proved  valuable  in  matching  simulator  performance 
to  measured  data. 

In  Figure  2.5  is  an  example  of  applying  the  simulator  to  a  specified  device  structure  with 
the  theoretical  doping  profile  is  shown  in  Figure  2.4.  The  device  was  a  C-band  MESFET 
fabricated  by  ITT.  The  device  has  a  gate  length  of  1  micron  and  a  gate  width  of  2.5  mm. 
The  specified  doping  profile  was  used  and  the  parasitics  and  electron  velocity-field  charac¬ 
teristic  were  adjusted  until  the  simulated  DC  I-V  characteristics  were  in  agreement  with  the 
measured  data.  As  indicated  in  Figure  2.5,  excellent  agreement  was  obtained.  In  order  to 
match  the  experimental  data  the  velocity-field  characteristic  for  GaAs  is  used  as  the  ‘fitting’ 
parameter.  The  low  field  mobility  can  be  taken  from  measured  data  or  varied  until  good 
agreement  between  the  measured  and  simulated  I-V  characteristic  in  the  linear  region  is 
obtained.  For  doping  densities  in  the  range  of  lO^^cm"®  the  mobility  is  usuaUy  in  the  range 
2000-4000  cm^/V -sec.  The  saturated  velocity  is  adjusted  until  the  model  predicts  the  correct 
open  channel  current  (la..)-  Typically,  for  a  one  micron  gate  length  device,  the  saturated 
velocity  is  in  the  range  of  1.4  —  1.6  x  10^  cm/sec.  The  saturated  velocity  scales  inversely 
with  gate  length,  as  expected. 

Simulations  of  the  class  A  RF  operation  of  the  device  at  5.5  GHz  produced  the  results 
shown  in  Figure  2.6.  As  indicated,  excellent  agreement  between  the  measured  and  simulated 
performance  is  obtained.  The  simulator  accurately  predicts  the  RF  output  power,  power- 
added  efficiency,  gain,  and  DC  drain  current,  over  a  range  of  input  RF  power  sufficient 
to  drive  the  device  better  than  6  dB  into  saturation.  Of  particular  note  is  the  agreement 
between  the  measured  and  simulated  PAE.  This  performance  function  is  typically  difficult 
to  match  using  commercially  available  equivalent  circuit  models. 

An  example  of  use  of  the  parameter  extraction  algorithm  is  shown  in  Figure  2.7.  This 
device  has  a  0.5  micron  gate  length  and  a  1  mm  gate  width  and  was  operated  class  A  at 
10  GHz.  The  specified  implant  profile  was  used  as  the  initial  doping  profile.  The  optimizer 
varied  the  doping  profile,  velocity-field  characteristic,  and  parasitics  until  good  agreement 
between  the  measured  and  simulated  performance  was  obtained.  Again,  excellent  agreement 
between  the  measured  and  simulated  data  was  obtained.  As  indicated  in  the  doping  profile 
plot,  the  free  electron  density  varies  from  the  implanted  donor  profile.  As  expected,  due 
to  diffusion  effects,  the  free  electron  distribution  is  reduced  near  the  implant  peak,  but  is 
increased  as  the  implant  ‘tails’  into  the  sub.strate.  Also  of  significance  are  the  variations  in 
the  drain  and  gate  currents  as  the  device  is  driven  into  saturation.  Accurate  modeling  of  the 
saturation  mechanisms  is  of  first  order  importance  to  obtaining  good  results.  It  has  been 
determined  that  saturation  in  MESFETs  occurs  due  to  forward  and  reverse  conduction  of 
the  gate  electrode.  This  will  be  discussed  in  more  detail  later  in  this  section. 

Two  other  ion-implanted  devices  operating  at  10  GHz  were  simulated  and  the  results 
obtained  are  shown  in  Figures  2.8  and  2.9.  Again,  excellent  agreement  between  the  measured 
and  simulated  data  is  obtained.  Devices  from  ITT,  Raytheon,  Texas  Instruments,  General 
Electric,  and  Hughes  have  been  simulated.  The  devices  had  gate  lengths  ranging  from 
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Figure  2.6:  ITT  MSAG  MESFET  -  Measured  and  simulated  performance  (Lg  =  1  fim,  Wg 
=  2.5  mm,  Va,  =  10  V,  la,  =  F  =  5.5  GHz) 
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0.25  micron  to  1  micron  and  gate  widths  ranging  from  0.5  mm  to  about  4  mm.  Both  ion- 
implanted  and  buried  channel  (i.e.,  low-high-low)  doping  profiles  were  simulated.  The  devices 
were  designed  for  operation  at  C  and  X-bands.  In  addition,  complete  integrated  amplifiers 
fabricated  by  ITT  were  simulated,  also  with  good  results.  Good  agreement  between  the 
simulated  and  measured  data  has  been  obtained  in  all  cases. 


2.4  Saturation  in  MESFETs 

The  excellent  agreement  between  the  measured  and  simulated  DC  and  RF  performance 
results  from  accurate  modeling  of  the  primary  saturation  mechanisms  in  MESFETs.  Investi¬ 
gation  of  saturation  in  these  devices  indicates  that  saturation  occurs  when  the  gate  electrode 
is  driven  into  reverse  and  forward  conduction.  This  is  indicated  in  Figure  2.10,  which  shows 
plots  of  the  dynamic  v-i  locus  superimposed  upon  the  DC  I-V  characteristics  for  a  typical  0.5 
micron  gate  length,  1  mm  gate  width  MESFET  operating  class  B.  The  dynamic  v-i  loci  are 
shown  for  operation  in  the  linear  region  and  when  the  device  is  driven  3  dB  into  saturation. 
When  operating  in  the  linear  region  the  dynamic  locus  is  essentisJly  elliptical,  except  when 
it  is  clipped  due  to  the  class  B  conditions.  As  saturation  is  achieved,  the  dynamic  terminal 
RF  voltage  exceeds  the  gate  breakdown  voltage  and  gate  current  flc'"s.  The  terminal  voltage 
decreases  until  it  gains  sufficient  magnitude  to  produce  forward  conduction  on  the  reverse 
portion  of  the  RF  cycle. 

When  operating  in  saturation  these  two  mechsmisms  compete.  That  is,  reverse  break¬ 
down  and  forward  conduction  occur  during  each  RF  cycle.  Which  mechanism  dominates 
is  determined  by  the  device  structure,  bias,  and  RF  tuning  conditions.  This  is  illustrated 
in  Figure  2.11,  which  shows  real-time  waveforms  at  the  gate  and  drain  terminals  for  a  0.5 
micron  gate  length  device  operating  class  A  at  10  GHz.  The  solid  lines  indicate  increasing 
levels  of  RF  drive,  varying  from  small-signal  linear  operation  to  about  3  dB  into  saturation. 
As  shown,  linear  operation  produces  sinusoidal  terminal  waveforms.  As  the  device  is  driven 
into  saturation,  significant  waveform  clipping  occurs.  The  waveform  distortion  is  primarily 
due  to  the  gate  conduction  characteristics,  as  indicated  in  the  two  plots  at  the  bottom  of 
the  figure.  In  this  particular  device,  forward  gate  conduction  is  dominant,  although  reverse 
breakdown  conduction  of  the  gate  electrode  is  also  clearly  seen. 

Reverse  breakdown  of  the  gate  electrode  is  dominant  in  determining  the  drain  voltage  at 
which  the  device  can  be  operated.  Typically,  for  optimum  performance  MESFETs  can  only 
be  biased  at  about  40-50%  of  the  drain-source  breakdown  voltage  due  to  the  saturation  mech¬ 
anism.  If  the  breakdown  voltage  is  low,  sufficient  drmn  bias  to  operate  the  device  efficiently 
can  not  be  applied  and  performance  suffers.  The  RF  output  power  that  can  be  obtzdned 
is  essentially  linearly  dependent  upon  drain  bias.  High  breakdown  voltage  is,  therefore,  of 
fundamental  importance  in  obtaining  maximum  RF  performance  from  MESFETs. 

Accurate  simulation  of  the  DC  and  RF  performance  of  MESFETs  requires  that  accurate 
models  for  gate  conduction  be  formulated.  Breakdown  in  MESFETs,  however,  is  not  well 
understood  and  an  accurate  model  that  predicts  observed  breakdown  phenomena  has  not 
previously  been  reported. 

Gate-drain  breakdown  must  be  accurately  simulated  if  good  agreement  between  mod¬ 
eled  and  experimental  data  is  to  be  obtained.  Unfortunately,  most  models  for  reverse  gate 
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Figure  2.8:  Measured  and  simulated  GaAs  MESFET  performance  (Lg  =  0.5  /xm,  Wg  =  1.25 
mm) 
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Figure  2.12:  I-V  characteristics  showing  drain-source  breakdown  for  ion-implanted  GaAs 
MESFET  (Lg  =  0.5  fim,  Wg  =  300  /im). 


breakdown  in  GaAs  MESFETs  do  not  yield  results  in  agreement  with  measured  data.  In 
particular,  the  models  do  not  correctly  predict  the  increasing  breakdown  voltage  versus  gate 
bias  shown  in  Figure  2.10  and  commonly  measured  in  experimental  data,  as  indicated  in 
Figure  2.12.  The  data  shown  in  Figure  2.12  were  measured  for  an  ITT  MSAG  MESFET. 
The  ability  to  simulate  this  behavior  is,  of  course,  fundamental  to  obtaining  accurate  sim¬ 
ulator  results.  We  have  proposed  a  mechanism  for  gate  breakdown  in  MESFETs  that  will 
yield  a  breakdown  model  suitable  for  incorporation  into  microwave  CAD  simulators  [17]. 
The  mechanism  has  been  verihed  with  experiment  [18].  A  simplified  and  qualitative  version 
of  this  model  is  partially  responsible  for  the  excellent  simulated  results  presented  here  by 
accurately  including  RF  saturation  due  to  gate  conduction.  The  model,  for  the  first  time,  cor¬ 
rectly  predicts  the  increasing  gate-source  breakdown  voltage  as  gate  bias  is  increased  towards 
pinch-off,  in  agreement  with  the  measured  data  shown  in  Figure  2.12.  The  gate  conduction 
model  allows  for  both  forward  and  reverse  breakdown  conduction  of  the  gate  electrode  and 
permits  calculation  of  conduction  during  the  RF  cycle,  as  shown  in  Figure  2.11. 
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Chapter  3 

Large-Signal  Analysis  and 
Optimization  of  Power  Amplifiers 


A  number  of  algorithms  were  added  to  the  simulator  to  calculate  large-signal  performance 
measures  applicable  to  MESFET  power  amplifiers.  These  large-signal  performance  measures 
were  then  used  as  the  basis  for  sensitivity  analysis,  nominal  performance  optimization,  and 
yield  optimization.  These  components  are  discussed  in  this  chapter  with  the  exception  of 
the  mathematical  optimization  algorithm,  which  is  the  topic  of  Chapter  4. 


3.1  Calculation  of  Performance  Measures 

Three  classes  of  performance  measures  are  available:  small-signal  gain,  measures  asso¬ 
ciated  with  gain  compression  levels,  smd  measures  associated  with  maximum  power-added 
efficiency.  Each  of  these  classes  has  two  subclasses.  The  first  subclass  consists  of  transducer 
quantities,  and  the  second  subclass  comprises  operational  quantities.  The  difference  between 
the  two  subclasses  will  be  developed  below. 

The  small-signal  transducer  gain,  Gf,  is  given  by: 


Where  Pg,dB  is  the  RF  power  in  dBm  available  from  the  generator  Vgen  shown  in  Figure  2.3 
and  Pi,dB  is  the  RF  power  in  dBm  delivered  to  the  load.  The  four  values  of  Pg,dB  “c  specified 
by  the  user,  and  they  must  be  in  the  small-signal  regime.  The  small-signal  transducer  gain 
has  no  other  associated  performance  measures. 

Similarly,  the  small-signsd  operation^  gain,  G„,  is  given  by: 

G.  =  T  Z  (PU  -  PtMB(Pi.6B))  ■ 

^  fc=l 

Where  Pi,dB  is  the  input  power  delivered  by  the  generator  to  the  circuit.  The  user  must 
supply  four  values  for  Pg,dB  that  are  all  in  the  small-signal  regime.  Given  these  values  of 
Pg,dBy  the  simulator  calculates  Pi,dB  and  Pi,dB-  The  smsdl-signal  operational  gain  has  no 
other  associated  preformance  measures. 
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The  transducer  gain  compression  measures  result  from  solution  of 

{Pl,dB{Pg,dB)  -  Pg,dB)  -  {Gt  -  Gc)  =  0  Pg,dB  €  R- 

Where  Gt  is  computed  by  the  formula  above  and  Gc  is  a  gain  compression  level  in  dB.  Gc  is 
user  selectable  as  1,  3,  or  6  dB.  The  Van  Wijngaarden-Dekker-Brent  root  finding  algorithm 
[19]  is  used  to  solve  the  above  equation.  Once  the  root  of  the  equation  is  found,  the  available 
performance  measures  are  Pg,dBi  Pi,dBi  and  the  gain-output  power  product  {2Pi,dB  —  Pg,dB) 
at  the  solution. 

The  operational  gain  compression  measures  result  from  solution  of 

{Pl,dB{Pg,dB)  -  Pi,dB{Pg,dB))  -  {Go  ~  Gc)  =  0  Pg,dB  ^  R- 

Where  Gc  is  the  selectable  gain  compression  level  and  Go  is  computed  by  the  previously 
given  formula.  Again,  the  Van  Wijngaarden-Dekker-Brent  root  finding  algorithm  [19]  is 
used.  Once  the  root  of  the  equation  is  found,  the  available  performance  measures  are 
PlABi  and  the  gain-output  power  product  (2Pj,<jb  —  P.y«)  at  the  solution. 

The  final  class  of  performance  measures  result  from  maximizing  expressions  for  the  power- 
added  efficiency.  The  maximum  transducer  power-added  efficiency  is  determined  as 

max  {PAEt}  =  max  I  f* 

where  Pi  is  the  RF  power  in  mW  delivered  to  the  load,  Pg  is  the  RF  power  in  mW  available 
from  the  generator,  and  Pdc  is  the  DC  power  in  mW  delivered  to  the  transistor.  Brent’s 
method  [19]  is  used  to  perform  this  1-D  optimization.  The  transducer  performance  measures 
are  the  maximum  PAEt  and  P„  Pj,  Pj/Pp,  and  PijPg  at  the  maximum. 

The  maximum  operational  power-added  efficiency  is  defined  as 

Brent’s  method  [19]  is  used  to  maximize  the  operational  power-added  efficiency  expression. 
The  operational  performance  measures  are  the  maximum  PAEo  and  Pj,  Pj,  Pi/Pi,  and  Pi  j  Pi 
at  the  maximum. 

To  calculate  the  performance  measures  the  user  must  specify  the  DC  bias  for  the  device. 
The  program  supports  eight  different  bias  modes  —  four  ways  of  setting  the  gate-source 
bias  voltage  and  two  ways  of  setting  the  drain-source  bias  voltage.  The  user’s  options  for 
specifying  the  gate-source  and  drain-source  bias  voltages  are  summarized  in  Table  3.1. 

While  calculating  a  performance  measure  a  number  of  situations  can  arise  where  the 
measure  become  invalid  for  some  reason.  These  situations  are  dealt  with  by  pen2dizing  the 
measure  under  consideration.  The  linear  value  of  the  measure  is  set  to  either  zero  or  1  x  10~*°. 
A  zero  value  is  used  if  the  linear  value  is  never  converted  into  a  logarithmic  representation. 
Otherwise,  the  value  of  1  x  10~*®  is  used.  The  circumstances  leading  to  penalization  of  the 
performance  measures  are  enumerated  below: 

1.  When  computing  gain  compression  or  maximum  power-added  efficiency  performance 
measures,  the  measure  is  penalized  if  the  gain  at  the  root  or  maximum  becomes  less 
than  1  dB,  or  some  other  user  specified  value. 
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Vg5  =  a 

User  specifies  a 

Vb5  =  bVg, 

User  specifies  b 

Idd{Vgs)  =  c 

User  specifies  c 

Ivd{Vgs)  =  dljDSS 

User  specifies  d 

11 

User  specifies  e 

^DS  =  f^dtbd 

User  specifies  / 

Table  3.1:  Bias  modes  available  when  calculating  performance  measures 

2.  The  performance  measure  is  penalized  when  the  device’s  pinch-off  voltage  fails  to  fall 
on  a  user  specifiable  interval,  which  defaults  to  -15  V  to  -1  V. 

3.  When  the  DC  bias  does  not  lie  in  the  valid  region  of  the  Vg,  —  V^  plane,  the  performance 
measure  is  penalized.  The  valid  region  is  given  by  the  following  relationships: 


V,, 

> 

AV^ 

default  A  =  0.95 

< 

BV^ 

default  B  =  0.1 

Vi, 

> 

CVi,u 

default  C  =  0.1 

Vi, 

< 

D  [mu(Vg,  —  Vpo)  +  Vi,u] 

default  D  =  1.0. 

Where  A,  B,  C,  and  D  are  user  controls  on  the  valid  region,  and  mu  is  slope  of  the 
line  in  the  Vj,  -  V.-  plane  which  marks  the  onset  of  gate-drain  breakdown. 

4.  Alternatively,  when  the  DC  bias  is  specified  as  an  invalid  —  Vd,  combination,  the 
performance  measure  is  penalized.  The  valid  /*i  —  V^,  region  is  given  by  the  following 
relationships: 

Idd  >  Eldt! 

Idd  <  FIdtt 

Vds  >  CVhm 

Vd,  <  D[mu{Vg,-Vg„)-\-V^M] 

Where  E  and  F  are  user  controls  on  the  valid  region. 

5.  The  device  and/or  its  embedding  circuit  are  supplying  RF  power  at  the  fundamental 
to  the  generator  impedance  (i.e.,  the  transistor  is  oscillating.) 

6.  While  evaluating  yield  estimates,  any  sample  devices  which  fall  outside  of  the  dis¬ 
turbance  hyperbox  (see  section  of  yield  estimation)  result  in  a  penalized  performance 
measure. 


default  E  =  0.05 
default  F  =  0.95 
default  C  =  0.1 
default  D  —  1.0. 


3.2  Large-Signal  Sensitivity  Analysis 


The  first  application  of  these  performance  measures  was  large-signal  sensitivity  analysis. 
The  sensitivity  is  given  by  the  following  formula: 


p  dm 
m  dp 
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where  m  is  the  performance  measure  and  p  is  the  parameter  under  study.  The  performance 
measures  available  are  all  those  listed  in  Section  3.1.  The  parameters  include  all  the  MESFET 
model  parameters  and  all  parameters  of  the  circuit  blocks  Zgen>  Z|o«d>  Zt,  and  Zfj,. 

The  sensitvity  is  calculated  numerically  within  TEFLON  by  evaluating  the  performance 
measure  at  a  sequence  of  user  specified  parameter  values.  This  sequence  of  data  points  is  then 
splined  with  a  Hermite  polynomial  spline.  The  spline  is  programmed  to  return  interpolates 
of  both  the  function  and  its  derivative.  The  derivative  interpolates  are  multilied  by  their 
corresponding  normalization  factors  Pifrrii  to  give  the  sensitivity. 

Sensitivity  analysis  is  useful  for  determining  which  device  or  circuit  parameter  fluctua¬ 
tions  cause  the  greatest  changes  in  critical  performance  measures.  This  knowledge  is  useful 
in  several  ways.  Once  a  nomined  set  of  device  and  circuit  parameters  is  selected,  sensitivity 
analysis  allows  assignment  of  parameter  tolerances  such  that  the  performance  measures  stay 
within  specified  bounds.  Another  application  of  sensitivity  analysis  is  determining  which 
device  or  circuit  parameters  to  use  during  optimization  runs. 

Examples  of  sensitivity  analyses  for  an  ion-implanted  device,  a  buried  channel  device,  and 
uniform  channel  device  are  given  in  Appendices  C,  D,  and  E,  respectively.  These  analyses 
were  performed  to  selected  the  optimization  variables  used  to  generate  the  results  presented 
in  Chapter  5. 


3.3  Nominal  Performance  Optimization 

Mathematical  optimization  of  a  device  design  (i.e.,  a  single  unperturbed  design)  with 
respect  to  one  of  the  above  performance  measures  is  nominal  performance  optimization.  At 
present,  nominal  performance  optimization  is  limited  to  a  single  measure  per  optimization 
run.  Multiple  performance  measures  could  be  used  if  a  satisfactory  method  of  combining 
them  into  a  single  number  was  selected  and  implemented. 

The  optimization  algorithm  employed  is  novel  and  was  developed  under  this  contract.  It 
is  a  gradient  based  algorithm  which  can  deal  successfully  with  objective  function  noise  and 
local  minimum. 

Nominal  performance  optimization  can  be  a  useful  design  tool,  but  often  the  resulting 
designs  exhibit  poor  manufacturing  yields.  These  poor  yields  result  from  the  optimizer 
placing  the  design  near  a  performance  cliff  [44,  42).  A  more  useful,  but  computationally  more 
expensive,  technique  is  yield  optimization,  which  was  the  main  effort  under  this  contract  and 
is  the  topic  of  the  next  section. 


3.4  Yield  Optimization 

Computer  aided  design  (CAD)  has  helped  to  improve  process  yields  and  the  average 
performance  of  monolithic  microwave  integrated  circuits  (MMICs)  through  the  application 
of  statistical  circuit  design  techniques.  A  principal  statistical  circuit  design  technique  is  yield 
optimization. 

The  yield  optimization  problem  can  be  formulated  two  ways: 

m«|T(a:)  =  j  p(v)(iv| 
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or 


mM|K(x)  =  j  p{v)4>{x  +  v)dv 

where  x  £  /T*,  p{v)  is  the  parameter  disturbance  probability  density  function,  and  Ra  is 
the  acceptability  region.  The  acceptance  function  <i>{x  -f-  v)  =  1  if  (x  +  v)  G  Ra-  Otherwise, 
^(x  +  v)  =  0.  The  first  formulation  requires  approximating  the  region  of  acceptability  and 
leads  to  region  of  acceptability  approaches.  The  second  formulation  is  usually  solved  with 
Monte  Carlo  techniques. 

A  number  of  different  authors  proposed  region  of  acceptability  approaches.  Scott  and 
Walker  [20]  and  Leung  and  Spence  [21]  pursued  the  regionalization  method.  Regionalization 
does  a  direct  search  on  the  space  to  determine  the  acceptabibty  region.  The  simplicial  ap¬ 
proach  was  used  by  Director  and  Hachtel  [22,  23,  24].  The  simplicial  approach  approximates 
the  region  of  acceptability  with  a  polyhedron  formed  by  points  on  the  boundary.  The  design 
center  is  then  the  center  of  a  hypersphere  which  is  contained  within  the  polyhedron.  Bandler 
and  Abdel-Malek  [25,  26,  27]  approximated  the  acceptability  region  with  hyperboxes  and 
linear  cuts.  The  yield  is  estimated  as  the  hypervolume  of  the  region  of  acceptability  divided 
by  the  hypervolume  of  the  tolerance  box.  The  multicircuit  approach  to  modeling  the  region 
of  acceptability  was  introduced  by  Bandler  et  al.  [28]  who  used  a  group  of  circuit  designs 
to  approximate  the  region  of  acceptability.  Bandler  and  Chen  [29]  then  used  generalized  Ij, 
centering  to  optimize  yield.  In  generalized  Ip  centering,  a  Ip  error  function  of  the  difference 
between  each  circuit’s  response  and  the  acceptable  region  boundary  is  minimized. 

Various  other  authors  sought  to  solve  the  yield  optimization  problem  with  Monte  Carlo 
techniques.  Soin  and  Spence  [30]  advocated  a  center  of  gravity  method.  N  Monte  Carlo 
samples  are  taken.  The  centers  of  gravity  of  the  passed  and  failed  points  are  determined. 
The  design  center  is  then  located  by  a  line  search  along  the  line  joining  the  centers  of  grav¬ 
ity.  Stochastic  optimization  was  used  by  Styblinski  and  Ruszczynski  [31]  and  Kjellstrom 
and  Taxen  [32].  In  this  approach,  small  sample  Monte  Carlo  jrield  estimates  are  stochas¬ 
tically  optimized.  Singhal  and  Pinel  [33]  introduced  parametric  sampling.  The  parametric 
sampling  technique  reuses  previous  samples  when  forming  a  Monte  Carlo  yield  estimate. 
The  control  variate  technique  [34]  for  reducing  the  Monte  Carlo  noise  in  the  yield  estimate 
was  used  by  Hocevar  et  al.  [35]  and  Soin  and  Rankin  [36,  37].  Biernacki  et  al.  [38]  used 
efRcient  quadratic  approximation.  This  technique  approximates  the  circuit  response  with 
a  multidimensional  quadratic  function  which  reduces  the  cost  of  the  Monte  Carlo  samples. 
The  approximation  is  generated  from  less  than  the  minimum  number  of  basis  points  by  a 
maximally  flat  interpolation  technique.  Ellipsoids  of  decreasing  volume  were  proposed  by 
Abdel-Malek  and  Hassan  [39].  This  method  approximates  the  region  of  acceptability  with 
an  ellipsoid  which  is  determined  by  decreasing  the  volume,  moving  the  center,  and  changing 
the  shape  of  a  starting  ellipsoid.  The  design  center  is  then  center  of  the  final  ellipsoid. 

In  addition  to  the  works  already  cited,  several  books  [40,  41,  42,  24,  43]  and  papers 
[29,  44]  review  this  area  of  research. 

The  previous  work  reported  on  yield  optimization  is  primarily  based  upon  the  use  of 
equivalent  circuit  models  for  active  devices.  This  dependence  ultimately  limits  present  gen¬ 
eration  CAD  since  equivalent  circuit  models  do  not  predict  operation  so  much  as  compactly 
represent  measurements  on  fabricated  devices.  Equivalent  circuit  elements  do  not  naturally 
capture  important  nonlinearities  and  interrelationships  of  the  physical  entities  they  repre- 
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sent.  Also,  determining  improved  device  design  parameters  from  equivalent  circuit  elements 
is  difficult  because  direct  correspondence  between  physical  parameters,  equivalent  circuit 
elements,  and  RF  performance  is  not  easily  established. 

In  this  work,  a  large-signal  GaAs  MESFET  simulator  for  yield  estimation  and  optimiza¬ 
tion  is  presented  that  does  not  rely  on  equivalent  circuit  techniques.  The  MESFET  model 
in  the  simulator  is  based  upon  device  physics  with  the  advantage  that  a  MIMIC  can  be 
simulated  from  process  data  all  the  way  through  to  RF  circuit  performance.  The  integrated 
simulator  allows  a  device  design  to  be  optimized  based  upon  a  desired  RF  performance  spec¬ 
ification.  The  small-  and  large-signal  power  amplifier  performance  measures  developed  in 
Section  3.1  are  available  for  the  yield  pass-faU  criterion.  The  yield  optimizer’s  variables  are 
physical  parameters  such  as  gate  dimensions,  channel  donor  distribution  specifications,  DC 
bias  voltages,  and  material  parameters.  A  standard  Monte  Carlo  method  is  used  to  predict 
the  process  yield  of  a  nominal  MESFET  design  based  upon  a  single  performance  measures’ 
variation.  The  Monte  Carlo  yield  estimate  is  then  optimized  using  a  quasi-Newton  method. 
The  quasi-Newton  method  is  deterministic  and  is  tolerant  of  the  inherent  noise  in  the  yield 
estimate. 

The  resulting  simulator  has  proved  to  be  accurate  under  both  class  A  and  B  operating 
conditions.  Since  simulations  can  be  performed  before  fabrication,  significant  time,  effort, 
and  expense  can  be  saved  in  the  development  of  advanced  MIMICs. 

Small  but  uncontrollable  disturbances  in  the  fabrication  process  result  in  devices  with 
geometries  and  doping  profiles  that  deviate  somewhat  from  nominal  values.  These  variations 
in  primary  process  parameters  appear  to  be  statistically  simpler  than  the  derivative  variations 
of  parameters  for  the  corresponding  equivalent  circriit.  For  example,  gate  width  and  length 
are  practically  independent  and  both  are  uncorrelated  with  ion  implant  dose  or  energy.  On 
the  other  hand,  variations  in  equivalent  circuit  parameter  values  such  as  g„,,  Rj,  Cg,,  Cdg, 
Cd»,  and  gd,  correlate  significantly  with  each  other  (45,  46,  47].  Moreover,  second  order 
correlations  (even  large  ones)  may  not  suffice  to  characterize  variations  in  equivalent  circuit 
parameters  [48]. 

Disturbances  in  physical  parameters  are  easily  characterized  by  a  second  order  statistical 
model,  a  multivariate  Gaussian.  A  multivariate  Gaussian  is  specified  by  a  mean  vector  and 
a  covariance  matrix.  The  mean  vector  is  simply  the  nominal  device  design.  The  covariance 
matrix  models  the  variances  and  covariances  between  physical  parameters. 

Given  a  Gaussian  model  for  the  parameter  disturbance,  the  yield  at  some  nominad  design, 
X,  is  estimated  by  the  following  Monte  Carlo  algorithm: 

1.  Input  mean  or  nominal  design  x, 

2.  Describe  disturbances  with  multivariate  probability  function  p{v), 

3.  Draw  a  set  of  N  disturbances  Vi  from  p{v),  and 

4.  Approximate  the  yield  as 

t 

The  disturbances  Vi  are  generated  as  follows.  Let  C  be  the  (positive  definite)  covariance 
of  the  parameters.  Solve  the  eigenvalue  problem  C  =  UWU^  for  the  matrix  of  eigenvectors  U 
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and  the  diagonal  matrix  of  (positive)  eigenvalues  W.  From  W  construct  the  diagonal  matrix 
w  with  entries  that  are  the  square  roots  of  the  entries  of  W.  Then  a  vector  Zi  of  independent, 
identically  distributed  Gaussian  numbers  of  zero  mean  and  unit  variance  is  constructed  by 
conventional  techniques  [19].  The  uncorrelated  disturbance  vector  Zi  is  converted  to  a  fully 
correlated  disturbance  vector  Vi  with  the  formula  Vi  =  Uwzi. 

The  acceptance  function  is  evaluated  at  each  vaJue  of  ®  by  comparing  the  per¬ 
formance  measure  under  consideration  with  a  user  specified  threshold,  (a;  -|-  Vi  is  shown  to 
have  the  correct  statistics  in  Appendix  B.)  If  the  performance  measure  is  greater  than  the 
threshold  value,  <f>  is  set  equal  to  one.  Otherwise,  <j)  is  set  equal  to  zero.  The  yield  estimate 
is  the  then  mean  value  of  the  acceptance  function  taken  over  the  N  samples. 

In  a  Monte  Carlo  yield  calculation,  many  devices  are  evaluated  in  the  vicinity  of  a  nominal 
design.  In  regions  where  the  performance  measure  is  relatively  smooth,  it  is  possible  to  use 
nearby  points  without  incurring  unacceptable  estimation  errors.  The  batch  size  of  a  Monte 
Carlo  calculation  sets  a  limit  on  expected  precision  which  can  be  related  to  derivatives  of 
performance  measure  to  determine  when  it  is  appropriate  to  replace  a  numerically  intensive 
performance  measure  calculation  with  a  previously  simulated  result. 

To  reduce  run  times  we  have  incorporated  a  simple  data  management  system  which 
maintains  a  binary  tree  database  of  previous  simulations.  When  a  new  design  is  to  be 
evaluated,  the  tree  is  searched  for  a  similar  design.  If  a  similar  design  is  close  enough  to  the 
desired  design,  the  stored  performance  measure  is  returned.  Otherwise,  the  device  model  is 
called,  and  the  simulation  result  is  recorded  in  the  tree. 

This  memory  system  will  reduce  the  number  of  calls  by  varying  amounts  dependent 
upon  the  particular  device  design,  the  degree  of  nonlinearity,  the  number  of  disturbance 
parameters,  and  the  extent  of  previous  data.  For  the  jrield  optimization  results  presented  in 
Chapter  5,  the  number  of  calls  was  reduced  by  16%. 
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Chapter  4 


Mathematical  Optimization 
Algorithm 


This  chapter  is  a  chronological  history  of  the  TEFLON  optimization  group’s  work  from 
June  1989  to  June  1990.  The  optimization  group’s  main  goal  was  to  assist  other  members 
of  the  TEFLON  team  with  problems  in  optimization.  The  group  also  helped  with  other 
mathematical  problems  that  the  TEFLON  team  encountered. 


4.1  Chronological  Program  Development 

The  first  task  for  the  optimization  group  was  to  obtain  and  test  various  optimization 
codes  from  the  netlib.  Two  codes  were  obtained,  LMDIF  from  Minpack  and  NL2SOL  from 
TOMS.  LMDIF,  a  least  squares  code,  proved  to  be  more  portable  and  was  more  successful  in 
solving  the  test  problems  we  ran  them  on.  LMDIF  used  the  Levenberg-Marquardt  algorithm 
with  a  forward  difference  approximation  to  the  Jacobian. 

During  the  summer  and  fall  of  1989  TEFLON  was  obtaining  poor  values  for  some  cal¬ 
culations,  particularly  near  the  end  points  of  intervals.  It  was  determined  that  part  of  the 
problem  was  due  to  the  interpolator  being  used  in  TEFLON.  We  wrote  a  cubic  piecewise 
Hermite  polynomial  interpolator  to  replace  the  existing  interpolator.  This  interpolator  uses 
a  five-point  algorithm  to  calculate  the  approximation  to  the  derivatives  at  the  interpolation 
points  for  the  desired  function.  The  interpolator  was  tested  against  several  codes  obtained 
from  the  netlib,  and  the  existing  interpolator  in  TEFLON,  which  was  a  code  from  IMSL. 
In  aU  the  problems  tested  our  interpolator  proved  to  be  superior.  In  particular,  it  obtained 
much  better  values  for  the  derivative  of  the  functions  to  be  interpolated  near  the  end  points 
of  the  interA^. 

In  the  spring  of  1990  the  group  began  working  on  device  optimization.  The  goal  in 
device  optimization  is  to  optimize  a  single  device  performance  function  such  as  input  power, 
output  power,  or  the  gain,  depending  on  from  two  to  thirty  parameters.  These  parameters 
include  gate  length,  gate  width,  peak  doping  density  in  the  channel,  etc.  These  device 
configuration  variables  must  fall  within  some  specified  bounds.  So  the  problem  becomes 
minimize  f  :  IV'  —*  R  ,  f  =  /(ii, •  •  •  ,x„)  under  the  constraints  that  li  <  Xi  <  /ti  , 
i  =  l,--,n. 
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Originally  we  were  attempting  to  use  least  squares  software  to  solve  this  problem.  To  use 
traditional  least  squares  techniques  we  attempted  to  solve  the  problem,  minimize  (C-f)  where 
C  >  maxf .  maxf  is  the  maximum  value  obtained  by  the  performance  function  /  on  the 
set  of  points  that  satisfy  the  above  inequality  constraints  on  the  parameters.  Unfortunately, 
the  code  the  group  was  using  (LMDIF,  from  Minpack)  would  not  converge  to  a  solution  if 
C  was  greater  than  maxf  by  more  than  a  relatively  small  amount.  Because  maxf  varied 
widely  as  a  function  of  the  type  and  number  of  the  parameters  used,  one  needed  to  know  the 
size  of  maxf  for  the  problem  one  was  working  on  to  a  very  close  approsdmation.  In  general 
this  was  not  the  case.  We  felt  that  all  least  squares  software  would  have  this  problem. 

We  next  tried  a  quasi-Newton  code  for  unconstriuned  optimization  to  work  the  device 
optimization  problem.  The  code  we  tried  used  the  BFGS  approximation  to  the  Hessian: 


Hk+i  =  Hk  + 


ykvl 

yl»k 


HkSkSkHk 

slHkSk 


where: 


=  ®fc+i  —  ®fc 


®fc+i  =  +  Pfc 

yk  =  V/(xfc+,)  -  V/(ifc). 

If  the  approximate  Hessian  is  not  positive  definite,  then  Hk+i  =  /• 

To  take  the  constraints  on  the  parameters  into  account,  the  algorithm  took  a  large  number 
of  aggressive  cutbacks  in  the  line-search.  Then,  if  the  algorithm  stepped  outside  of  the 
hyperbox  defined  by  the  constraints,  it  was  hoped  that  it  would  be  able  to  quickly  backtrack 
into  the  allowable  area.  However,  when  the  code  stepped  outside  of  the  hyperbox,  numerical 
instability  was  encountered  with  the  TEFLON  software  package. 

We  then  tried  working  the  problem  using  a  projected  gradient  code.  This  code  is  a  steep¬ 
est  descent  code  as  long  as  the  steps  stay  inside  of  the  hyperbox  defined  by  the  constraints 
on  the  parameters.  If  the  code  steps  outside  of  the  hyperbox  the  step  is  the  projection  of 
the  steepest  descent  step  onto  the  surface  of  the  hyperbox.  The  projection  of  the  steepest 
descent  step  onto  the  hyperbox  is  accomplished  by  the  following  algorithm: 


If  Xi  —  /3"‘V/(x*)i  <  li,  then 

X?  =  ti. 

If  X*  —  /?”‘V/(x*),  >  fii,  then 
=  fii- 
Otherwise 

=  xf  -  P”^Vf{x%. 
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Where: 

1.  if  is  the  component  of  the  k‘^  point  generated  by  the  algorithm. 

2.  li  and  /x,-  are  the  lower  and  upper  bounds  on  the  i**'  parameter,  respectively. 

3.  V/(a;*')i  is  the  i‘^  component  of  the  gradient  of  f  at  the  k‘*'  point. 

4.  yd  is  a  user-supplied  constant. 

5.  m  is  the  first  nonnegative  integer  such  that 

/(i^+i)  <  /(**)  +  -  I*) 

is  satisfied,  o’  is  a  user-defined  constant. 

6.  /  is  scaled  to  be  approximately  unity. 

Although  this  code  was  more  successful  then  amy  code  previously  tried,  it  still  became 
trapped  in  local  minima.  Thus,  it  was  unable  to  obtain  the  global  minima. 

In  an  attempt  to  understand  the  noise  in  the  problem,  we  tabulated  values  of  the  per¬ 
formance  measure,  small  signal  gain,  for  points  on  a  rectangular  grid  covering  the  allowable 
parameter  space  for  several  pairs  of  parameters.  The  resulting  tables  indicated  two  levels  of 
structure.  The  larger  level  was  a  single  basin  with  its  minimum  lying  somewhere  within  the 
hyperbox.  The  finer  level  of  structure  imposed  a  rough  surface  on  the  basin.  This  surface 
had  many  local  minima.  It  was  these  local  minima  that  were  trapping  the  projected  gradient 
algorithm. 

In  June  1990,  the  optimization  group  implemented  the  first  version  of  the  current  code. 
This  new  implementation  of  the  projected  gradient  algorithm  used  a  monotonically  decreas¬ 
ing  sequence  of  steps  for  the  finite  difference  approximation  to  the  gradient.  The  initial  step 
was  taken  to  be  half  the  length  of  the  shortest  side  of  the  hyperbox.  With  this  size  of  step 
the  descent  direction  given  by  the  approximate  gradient  was  more  affected  by  the  global 
structure  of  the  performance  measure  than  by  the  local  structure  near  the  current  iteration 
point.  This  way  the  algorithm  was  able  to  avoid  the  many  local  minima  caused  by  the  rough 
surface  and  move  the  iteration  points  toward  one  of  the  minima  of  the  global  structure. 

A  finite  difference  step  size  was  rejected  under  two  conditions: 

1-  If  lIV/lloc  <  iol.  Note  that  V/  is  the  finite  difference  approximation  to  the  gradient. 
tol  is  some  estimate  of  the  level  of  noise  in  the  performance  measure. 

2.  The  backtracking  line-search  fails  to  find  a  point  on  the  line  segment  between  and 
I*  —  V/(a:*)  that  satisfies  the  Armijo  condition. 

If  either  of  the  above  two  conditions  held,  the  finite  difference  step  size  was  cut  by  two  and 
the  algorithm  continued,  starting  from  the  last  point  generated  by  the  algorithm  with  the 
previous  finite  difference  step  size.  Decreasing  the  step  size  in  this  manner  allowed  the  al¬ 
gorithm  to  work  its  way  down  into  one  of  the  basins  of  the  global  structure.  The  algorithm 
terminated  when  the  step  became  so  small  that  the  finite  difference  gradient  was  only  mea¬ 
suring  changes  in  the  rough  surface.  This  algorithm  proved  successful  in  several  test  problems 
with  similar  structure  to  the  device  optimization  problems  we  were  currently  working  on. 
More  conventional  optimization  techniques  were  unable  to  find  the  global  minima  of  these 
test  problems.  The  new  algorithm  was  adso  tested  on  several  device  optimization  problems 
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using  the  TEFLON  software  package.  In  each  of  these  problems,  the  algorithm  was  started 
from  a  number  of  different  initial  points.  Given  a  problem,  the  algorithm  converged  to  the 
same  global  minimum  from  a  variety  of  starting  points.  The  algorithm  was  generally  able 
to  satisfy  ||V/|joo  <  tol  at  this  point,  for  the  finite  difference  step  size  not  too  small.  The 
algorithm  was  tested  against  a  simulate  annealing  code  on  two  dimensional  problems.  Gen¬ 
erally,  the  algorithm  achieved  larger  values  for  the  performance  measure  and  it  converged 
within  far  fewer  function  evaluations. 

In  July  the  single  finite  difference  step  was  replaced  with  a  vector  of  finite  difference  step 
sizes.  Each  element  in  this  vector  corresponded  to  one  of  the  parameters  in  the  problem. 
The  initial  finite  difference  step  for  each  component  of  the  gradient  was  set  to  half  the  length 
of  the  side  of  the  hyperbox  corresponding  to  that  parameter.  A  vector  of  finite  difference 
steps  was  rejected  under  the  same  conditions  as  a  finite  difference  step  in  the  first  code.  The 
vector  of  steps  was  then  multiplied  by  0.5  and  the  algorithm  started  again. 

Three  modifications  were  made  to  the  code  in  the  month  of  August. 

•  If  II V/ II  <  1 1  he  1 1  set  II  V/ll  =  II  he  II .  Where  he  is  the  current  finite  difference  step  size  and 
II  •  II  indicates  the  Euclidean  norm.  We  tested  this  modification  on  two  problems.  The 
two  problems  used  for  the  experiment  both  used  small  signal  gain  for  the  performance 
measure,  and  the  following  parameters: 

-  Nmax,  1  X  10^®  <  Nmax  <  1  x  10^*. 

-  Range  ,  0.4  <  Range  <  0.12. 

with  the  starting  point  Xo  =  (Nmax,  Range)  =  (9  x  10^^,  0.10).  However,  the  applied 
voltage  was  calculated  differently  in  the  two  problems.  The  two  methods  were: 

-  Method  1:  V/  unrestricted. 

-  Method  2:  ||V/||  >  ||h||  where  h  is  the  vector  of  finite  difference  step  sizes. 

The  following  table  was  obtained: 


method  1 

method  2 

PM 

NFE 

PM 

NFE 

A 

1.282 

66 

1.409 

63 

B 

1.879 

33 

1.879 

24 

In  the  table  PM  stands  for  performance  measure,  and  NFE  stands  for  number  of 
function  evaluations.  From  the  above  table  it  appears  that  method  2  is  superior.  From 
plots  of  the  surfaces  of  the  problems  defined  above,  the  values  obtained  in  the  table 
were  seen  to  be  global  maxima.  In  another  problem  not  listed  in  the  table  method  2 
was  able  to  converge  to  a  solution  while  method  1  was  not. 

•  Default  centered  differences  was  replaced  with  default  forward  differences.  Four  ex¬ 
periments  were  run  each  with  small  signal  gain  as  the  performance  measure.  The 
experiments  are  summerized  below: 
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1  Run  1 

Parameters 

Bounds 

Starting  Point 

IIB 

Gate  length,  Sigma 

[0.2, 1.0],  [0.1, 0.2] 

(0.9,0.19) 

B 

Sigma,  Range 

C 

Nmsix,  Range 

1  X  101«,1  X  10^«],  [0.04,0.12] 

(9  X  10l^0.10) 

D 

Nmax,  Gate  Width 

[1  X  10^«,1  X  10^®],  [500.0, 2000.0] 

The  two  methods  used  were: 

—  Method  1:  V/  obtained  using  forward  differences. 

—  Method  2:  V/  obtained  using  centered  differences. 

The  following  table  was  obtained. 


method  1 

method  2 

PM 

NFE 

PM 

NFE 

A 

31.31 

60 

31.06 

53 

B 

14.43 

58 

14.68 

88 

c 

14.09 

63 

14.10 

74 

D 

15.81 

62 

17.13 

76 

In  the  table  PM  stands  for  performance  measure,  and  NFE  stands  for  number  of 
function  evaluations.  From  the  above  table  it  appears  that  using  centered  differences 
gives  a  more  robust  algorithm  than  using  forward  differences. 

•  The  projected  gradient  step  was  replaced  with  a  projected  BFGS  step.  Three  methods 
were  tested  against  each  other. 

—  Method  1:  Projected  gradient. 

—  Method  2:  Projected  BFGS,  where  the  approximate  Hessian  is  set  to  the  identity 
every  time  a  scale  is  rejected. 

-  Method  3:  Projected  BFGS,  where  the  approximate  Hessian  is  updated  through¬ 
out  the  run. 

The  Hessian  was  set  to  the  identity  whenever  it  was  not  positive  definite  in  both 
projected  BFGS  methods.  The  two  problems  used  were: 


Run 

Parameters 

Bounds 

Starting  Point 

A 

Sigma,  Range 

[0.1, 0.2],  [0.04,0.12] 

(0.19,0.11) 

B 

Nmax,  Range 

[1  X  10^®,!  X  10'«],  [0.04, 0.12] 

"(9TTF^;0)" 

The  following  table  summarizes  the  results. 


method  1 

method  2 

method  3 

NFE 

PM 

NFE 

PM 

D 

88 

14.68 

81 

14.63 

55 

12.22 

B 

63 

14.09 

65 

13.41 

63 

14.15 
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The  above  table  does  not  give  conclusive  evidence  as  to  which  is  the  best  method. 
More  experimentation  is  necessary. 

During  the  fall  of  1990  we  began  working  problems  of  interest  to  other  members  of  the 
TEFLON  team.  The  first  problem  we  worked  was  to  maximize  the  gain  for  ITT  169  wafer 
MESFETs  being  operated  in  the  class  B  mode.  This  was  to  be  done  by  varying  the  channel 
doping  profile.  Several  tests  were  run  on  models  of  the  above  mentioned  devices.  The  code 
was  able  to  increase  the  gain  by  10%  to  15%  over  the  gain  obtained  by  manually  adjusting 
the  parameters.  However,  in  doing  so  the  code  drove  the  device  model  out  of  the  class  B 
mode  and  into  the  class  A  mode.  The  TEFLON  source  code  was  altered  so  that  the  user 
could  specify  what  mode  the  model  should  operate  in.  With  these  changes  in  place,  the 
code  was  able  to  increase  the  performance  measure,  maximum  power-added  efficiency  (Max. 
PAE)  depending  on  the  three  channel  doping  profile  parameters  nmax,  range,  and  sigma,  by 
14%  over  a  device  obtained  by  manually  adjusting  the  parameters.  This  improvement  was 
obtained  by  using  the  parameters  of  the  manually  optimized  device  as  a  starting  point. 

However,  when  using  other  starting  points  the  code  found  other  local  minima  with  im¬ 
possibly  large  values  of  Max.  PAE.  After  looking  at  a  number  of  graphs  of  the  performance 
measure  on  two  dimensional  subspaces  of  the  parameter  hyperbox,  it  was  decided  that  the 
impossible  performance  measures  were  caused  by  impossibly  large  gate  to  source  voltages. 
A  penalty  function  was  added  to  the  performance  measure  that  helped  keep  the  algorithm 
out  of  regions  of  the  hyperbox  where  these  impossibly  large  gate  to  source  voltages  occurred. 

At  this  time  the  group  began  working  on  several  more  problems: 

1.  Maximize  maximum  power-added  efficiency  for  ITT  169  class  B  MESFETs  by  adjusting 
the  matching  circuit  parameters:  the  real  and  imaginary  parts  of  the  load  and  generator 
impedances  at  the  fundamental,  second,  and  third  harmonics. 

2.  Maximize  small  signal  gain  for  ITT  169  MESFETs  by  adjusting  the  matching  circuit 
parameters. 

3.  Maximize  maximum  power-added  efficiency  for  ITT  169  class  B  MESFETs  by  adjusting 
the  doping  profile  parameters  and  the  matching  circuit  parameters. 

The  first  problem  was  of  particular  interest  at  the  time.  The  algorithm  was  able  to  obtain 
the  same  high  performance  measure  (Max.  PAE>  80%)  from  many  different  starting  places 
spread  around  the  hyperbox.  All  of  the  final  iteration  points  were  within  a  few  ohms  of  one 
another.  The  algorithm  was  able  to  obtain  this  high  performance  measure  within  several 
hundred  function  evaluations.  This  compared  very  favorably  to  the  simulated  annealing  code 
also  being  used  in  TEFLON. 

During  the  month  of  November  several  changes  were  made  to  the  code  for  theoretical 
reasons.  These  changes  were: 

•  Scaling  -  the  hyperbox  was  mapped  to  the  unit  cube. 

•  After  going  through  all  scales  the  algorithm  started  over  again  at  the  first  scale,  con¬ 
tinuing  until  no  progress  was  made  at  any  scale. 
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•  The  method  for  determining  convergence  at  a  particular  scede  was  changed  from 
liv/ll  ^  tol  (where  iol  is  an  estimate  of  the  noise  in  the  performance  measure  )  to 

iiv/ii  <  S 

The  code  with  these  three  changes  was  never  tested  against  the  previous  version  of  the 
code.  These  changes  were  made  for  theoretical  reasons.  However  a  slight  improvement  in 
the  performance  of  the  algorithm  seemed  to  be  obtained. 

During  the  month  of  December  various  stopping  criteria  were  tested.  These  criteria  were 
equivalent  to  “convergence”  at  the  current  scale.  All  tests  for  convergence  were  made  after 
the  gradient  at  the  current  point  was  set  to  at  least  as  large  as  the  current  scale.  Four 
different  criteria  were  tested: 

•  criteria  1:  ||V/||  <  where  he  is  the  current  finite  difference  step  and  tol  is  an 
estimate  of  the  noise  in  the  function. 

•  criteria  2:  ||V/||  <  O.lhe. 

•  criteria  3:  jj®  —  (®  —  V/)'^|)  <  where  A  indicates  the  projection  onto  the  hyperbox. 

•  criteria  4:  ||®  —  (®  —  V/)^||  <  O.lhc. 

Six  experiments  were  used  for  the  tests.  AU  of  the  experiments  were  done  on  a  ITT  169  class 
A  MESFET  and  had  small  signal  gain  as  the  performance  measure.  The  following  table 
summerizes  these  six  problems. 


Run 

Parameters 

Bounds 

Starting  Point 

A 

Nmax,  Sigma 

[1  X  10*^6  X  10'^],  [0.05, 0.30] 

B 

Nmax,  Sigma 

[1  X  10‘^6  X  10»n,  [0.05, 0.30] 

graaaniMiiEiiii 

C 

Nmax,  Sigma 

[1  X  10l^6  X  10^^],  [0.05, 0.30] 

(4.6  X  10^^0.109) 

D 

Nmax,  Range 

[1  X  10^^ 6  X  10^^],  [0.001, 0.20] 

"(IXxlo^’ToioT^ 

E 

Nmax,  Range 

[1  X  10^^6  X  10^7],  [0.001,0.20] 

(3.6  X  10^^0.179) 

mm 

Nmax,  Range 

The  results  of  these  tests  are  summarized  in  the  following  table: 


criteria  1 

criteria  2 

criteria  3 

criteria  4 

NFE 

PM 

NFE 

PM 

NFE 

PM 

NFE 

PM 

A 

77 

13.171 

77 

13.171 

77 

13.171 

77 

13.171 

B 

76 

13.206 

98 

13.207 

76 

13.206 

98 

13.207 

C 

61 

13.204 

82 

13.208 

61 

13.204 

82 

13.208 

D 

72 

13.139 

88 

13.140 

72 

13.139 

93 

13.140 

E 

71 

13.140 

no 

13.140 

71 

13.140 

142 

13.140 

F 

54 

12.006 

DNC 

DNC 

44 

12.006 

DNC 

DNC 

DNC  means  the  algorithm  did  not  finish  within  ten  minutes  of  CPU  time. 

As  can  be  seen  from  the  above  table,  rejecting  the  finite  difference  step  size  when 

IN  -  (*  -  v/)-'!!  < 
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appears  to  be  the  best  method. 

In  January  1991  the  code  was  modified  further.  Several  changes  were  tried.  The  first  was 
moving  the  check  for  convergence  from  after  the  gradient  is  extended  to  before  the  gradient 
is  extended.  The  second  wa.s  changing  the  projected  gradient  code  to  a  projected  SRI  code. 
The  formula  for  calculating  the  SRI  approximation  to  the  Hessian  is  given  below: 


Sk+i 


Where  5*  is  the  previous  SRI  approximation  to  the  Hessian, 


a*  =z  —  a;* 


and  r*  =  where  y*  =  V/(®*+^)  -  V/(®*). 

The  performance  measure  used  in  these  experiments  was  small  signal  gain.  Eight  experiments 
were  run.  They  are  summarized  in  the  following  table: 


Run 

Problem 

Bounds 

Starting  Point 

A 

Nmax,  Range 

[1  X  101^6  X  10‘^],  [0.001, 0.20] 

B 

Nmax,  Range 

C 

Nmax,  Range 

[1  X  10*’^,  6  X  10‘^],  [0.001, 0.20] 

D 

Nmax,  Sigma 

[iTiF^;?iriFn7[o:o5roo] 

(1.6  X  10”,  0.109) 

E 

Nmax,  Sigma 

(3.6  X  10”,  0.199) 

F 

Nmax,  Sigma 

[1  X  10‘S6x  10”],  [0.05, 0.30] 

jT6no^'’70o9j 

G 

Nmax,  Sigma,  Range 

■fnnriUiiimiuiHOkiioiiijFiiiAfiiiidtiMiji 

(3.6  X  10*^0.079,0.199) 

H 

Nmax,  Sigma,  Range 

[1  X  10”,  6  X  10‘'J,  [0.001, 0.20],  [0.05, 0.30] 

(4.0  X  10”,  0.129, 0.210) 

The  two  methods  used  were: 

•  Method  1:  projected  gradient. 

•  Method  1:  projected  SRI. 

The  results  of  these  experiments  are  summerized  in  the  following  table. 


method  1 

method  2 

PM 

NFE 

PM 

NFE 

D 

13.139 

72 

B 

13.140 

71 

13.140 

71 

C 

12.006 

44 

12.006 

44 

D 

13.171 

77 

13.207 

75 

E 

13.206 

76 

13.205 

mm\ 

H 

13.204 

61 

13.204 

56 

G 

13.781 

97 

13.810 

86 

H 

12.873 

218 

13.004 

188 

In  experiment  C  nether  algorithm  was  able  to  move  from  the  initial  point.  In  experiment  H 
the  number  of  cutbacks  was  set  to  10  instead  of  4  as  in  the  other  experiments.  As  can  be 
seen  from  the  above  table  SRI  did  as  well  or  better  on  every  problem. 
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During  the  month  of  March  a  cubic  model  of  the  performance  measure  in  the  quasi- 
Newton  direction  was  added  to  the  line-search.  See  section  6.3.2  of  Numerical  Methods  for 
Unconstrained  Optimization  and  Nonlinear  Equations,  by  Dennis  and  Schnabel  for  details 
of  the  implementation.  This  model  is  used  to  calculate  the  optimum  value  for  the  cutback 
factor.  One  deviation  was  made  from  the  implementation  as  given  by  Dennis  and  Schnabel. 
The  line-search  was  terminated  after  10  cutbacks  and  not  when  the  size  of  the  reduced  step 
became  less  than  some  specified  value.  This  method  of  doing  the  line-search  was  compared 
to  the  linear  cutback  strategy  of  using  cutbacks  of  |  that  wais  currently  being  used  in  the 
algorithm.  The  performance  measure  used  in  these  experiments  was  small  signal  gain.  Eight 
experiments  were  run.  They  are  summarized  in  the  following  table: 


Run 

Problem 

Bounds 

Starting  Point 

A 

Nmax,  Range 

[1  X  101^6  X  10^^,  [0.001, 0.20] 

(1.6  X  lo"^  0.079) 

B 

Nmax,  Range 

[1  X  101^,6  X  10*^],  [0.001, 0.20] 

Nmax,  Range 

[TTTF'TeTlF’tTlOoIToISo] 

D 

Nmax,  Sigma 

E 

Nmax,  Sigma 

(3.6  X  10‘', 0.199) 

F 

Nmax,  Sigma 

[1  X  101^6  X  lOi’^],  [0.05, 0.30] 

(4.6  X  10”,  0.109) 

Nmax,  Sigma,  Range 

[1  X  10‘',6  X  10‘'J,  [0.05,0.30],  [0.001,0.20] 

“(aTxTF'TolisoToTo^ 

H 

Nmax,  Sigma,  Range 

“(3T>rTF7o^20M^oT)“ 

In  experiments  C,  G,  and  H  the  minimum  scale  was  set  to  0.00005  and  the  number  of 
cutbacks  used  was  10.  In  the  remaining  experiments  the  minimum  scale  was  set  to  0.005, 
and  the  number  of  cutbacks  used  was  4.  In  all  cases  the  scaling  factor  for  the  performance 
measure  was  The  two  methods  are  listed  as; 

•  Method  1:  Projected  SRI  with  cubic  line-search. 

•  Method  2:  Projected  SRI  with  linear  line-search. 

The  resvilts  of  these  experiments  are  summerized  in  the  following  table. 


method  1 

method  2 

PM 

NFE 

PM 

NFE 

a 

13.138 

59 

13.140 

66 

B 

13.140 

80 

13.140 

71 

c 

13.140 

235 

13.140 

259 

D 

13.207 

kb 

13.207 

75 

E 

13.207 

73 

13.205 

75 

F 

13.207 

59 

13.204 

56 

G 

13.107 

201 

13.170 

246 

H 

13.327 

293 

13.327 

316 

As  can  be  seen  from  the  above  table  the  cubic  model  performed  better  in  6  out  of  the  8 
experiments.  It  also  did  better  in  every  case  where  a  large  number  of  function  evaluations 
were  done.  Because  of  these  results  the  cubic  model  will  be  the  line-search  procedure  used 
from  now  on. 

In  April  the  group  assisted  the  engineers  in  the  TEFLON  group  with  applying  the  algo¬ 
rithm  to  their  work.  This  involved  looking  at  the  output  from  runs  on  their  problems  and 
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advising  them  how  to  set  the  various  parameters,  such  as  the  scaling  factor  for  the  perfor¬ 
mance  measure,  the  minimal  size  for  the  finite  difference  steps  and  the  number  of  cutbacks 
to  be  used  in  the  line-search.  This  work  was  useful  not  only  because  it  helped  the  engineers 
accomplish  their  work,  but  it  also  gave  us  insights  into  the  performance  of  the  code. 

We  also  began  working  on  the  theoretical  aspects  of  the  algorithm.  This  work  involved 
analyzing  the  behavior  of  projected  gradient  algorithms  in  the  presence  of  noise  in  the  perfor¬ 
mance  measure  and  when  analytical  derivatives  are  not  available.  We  feel  that  this  work  will 
give  insight  into  the  workings  of  the  algorithm.  In  particular  it  is  hoped  that  this  analysis 
will  give  insight  into  when  to  terminate  the  iteration,  and  on  rates  of  convergence. 


4.2  Current  Algorithm 


A  description  of  the  current  algorithm  is  as  follows: 

1.  Given 


calculate 

If 

set 


f{x‘),  V/(x'),  and  y  =  ix‘-  V/(x'))\ 


ll*‘-y||  < 


minacal 

2h 


h  =  O.bh,  and  5'  =  / 


then  return  to  1.  (V)''  indicates  the  projection  of  V  onto  the  hyperbox,  and  minacal  is 
the  minimal  scale  used,  minacal  is  an  estimate  of  the  minimum  distance  between  points 
in  the  hyperbox  for  which  an  appreciable  difference  in  the  value  of  the  performance 
measure  can  be  detected  and  5'  is  the  SRI  matrix. 


2.  Otherwise  .'pda  e  5'  using  the  SRI  update  described  earlier  in  this  report.  Determine 
if  S‘  is  positive  definite.  If  it  is  not,  set  5'  =  /.  In  our  code  positive  definiteness 
of  the  matrix  is  checked  in  the  Linpack  code  DCHDC,  a  double  precision  Cholesky 
decomposition  routine.  Setting  =  /  if  the  approximate  Hessian  is  not  positive 
definite  guarantees  that  the  step  is  always  a  descent  direction. 


3.  Solve 

5'p  =  V/(x'). 


If 


IIpIU  <  h 


then  set 


p  = 


h 


The  matrix  equation  is  solved  by  use  of  the  Linpack  subroutines  DCHDC  and  DPOSL. 
First  DCHDC  does  a  Cholesky  decomposition  on  the  matrix  and  then  DPOSL  solves 
the  factored  system. 
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4.  Calculate  the  cutback  factor  a  for  the  step.  The  algorithm  for  calculating  a  will  be 
discussed  latter  in  this  report,  below  for  details.  If 


set 

then  return  to  1. 
5.  Calculate 

If 

then  set 


\\ap\\  <  h, 

h  =  0.5/1,  and  5'  =  / 

=  (x'  -  Qp)^. 

f(z*)  <  /(*')  -  10-<aV/V 


x‘  =  x'*' 


and  return  to  1.  Otherwise  return  to  4. 

The  algorithm  for  calculating  the  cutback  factor  for  the  line-search  goes  as  follows; 
If  no  previous  cutbacks  for  this  p  have  been  made,  set  a  =  1.  If 

X;  =  Ui  and  x‘  ^  Ui  or  xf  =  li  and  x'  ^  U 

for  any  i  set 

a  =  0.5a. 

Where  a  is  the  previous  cutback  factor  that  was  used  to  calculate  x+ .  If  q  is  the 
first  cutback  factor  for  this  p  such  that 

x^  ^  Ui  if  x'  ^  Ui,  and  xf  ^  li  if  x‘  ^ 

for  every  t,  then  a  should  be  the  unique  minimizer  of 

(/(x‘  -  ap)  -  /(x')  -  dV/(x')^p)a’  +  dV/(x')^pQ  -f  /(x'). 

Otherwise,  a  is  the  local  minimizer  to 

OQ^  -I-  ba^  +  V  f{x‘)^ pa  -f  /(x'). 

Where: 


a 

1 

■  1 

-1  1 
F 

f{x‘ 

-  ap)  -  /(*")  -  ^f{x‘)‘^pa 

b 

a  —  a 

d 

r  J 

/(x' 

-  ap)  -  /(*  )  -  '^/(®')^pa 

In  the  above  matrix  equation,  d  is  the  cutback  factor  previous  to  d.  The  local 
minimizer  to  this  cubic  polynomial  is 


a  = 


-6-f  yjb^  ~  ZaV f{x‘yp 


3a 
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The  code  as  described  above  has  been  successful  in  solving  a  number  of  problems  of 
interest  to  the  TEFLON  team.  These  problems  include: 

•  Device  optimization:  modifying  the  physical  parameters  of  the  device  so  as  to  obtain 
maximum  performance  for  a  given  performance  measure. 

•  Yield  optimization:  centering  the  design  of  a  device  within  the  hyperbox  defined  by 
the  physical  parameters  to  obtain  the  greatest  number  of  manufactured  devices  that 
perform  above  a  given  threshold  for  a  given  performance  measure. 

•  Parameter  extraction:  modifying  the  physical  parameters  so  as  to  match  the  perfor¬ 
mance  of  a  given  given  device  for  a  given  set  of  performance  measures. 
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Chapter  5 

Yield  Optimization  Experiments 


The  simulator  was  used  to  conduct  a  series  of  yield  optimization  experiments.  Three  types 
of  MESFET  devices  were  selected:  an  ion  implanted  device,  a  buried  channel  device,  and  a 
uniform  channel  device.  The  ion  implanted  MESFET  design  was  based  upon  a  commercial 
device  that  had  been  empiricsdly  optimized  by  a  standard  procedure  for  maximum  power- 
added  efficiency.  The  device  is  capable  of  good  RF  performance  and,  as  will  be  shown, 
the  simulator  was  not  able  to  alter  the  design  to  obtain  significant  design  improvements  for 
PAE.  The  buried  channel  MESFET  design  was  also  based  upon  an  industrial  device,  but 
the  device  in  this  case  was  a  prototype  and  was  not  optimum.  The  simulator  significantly 
modified  the  original  design  to  improve  performance.  The  starting  design  for  the  uniform 
channel  MESFET  was  determined  using  standard,  first  principle  design  techniques.  The 
initial  design  proved  not  to  be  optimum  and  was  also  significantly  modified  by  the  simulator 
to  obtain  improved  performance.  Starting  from  these  MESFET  designs,  we  then  sought  to 
optimize  the  yield  of  each  device  structure  relative  to  to  small-signal  transducer  gain,  output 
power  at  1  dB  gain  compression,  and  maximum  power-added  efficiency.  All  devices  were 
embedded  in  a  50  H  circuit  so  that  performance  variations  based  upon  device  design,  rather 
than  circuit  tuning  conditions,  could  be  investigated. 

The  experiments  were  conducted  by  specifying  the  experiment  variables  and  estimating 
the  disturbance  covariance  matrices.  We  selected  the  variables  by  performing  sensitivity 
analyses  on  each  device  structure.  The  results  of  these  analyses  are  given  in  Appendices  C, 
D,  and  E.  The  variables  to  which  the  performance  measures  are  most  sensitive  were  included. 
For  all  the  device  types,  the  performance  measures  were  most  sensitive  to  the  variables  which 
specify  the  gate  geometry,  the  channel  donor  distribution,  and  the  DC  bias.  The  variables  for 
each  device  type  are  listed  in  Table  5.1.  For  these  experiments  the  gate  length  is  perturbed 
about  the  nominal  manufacturable  value  which  held  fixed  and  is  not  subject  to  optimization. 
Previous  experience  indicates  that  the  optimizer  always  drives  gate  length  to  its  minimum 
permitted  value. 

We  sissumed  that  each  device’s  variables  are  statistically  independent.  The  covariance 
matrices  are  therefore  diagonal.  The  diagonal  elements  are  the  squares  of  estimated  standard 
deviations  for  each  variable.  In  all  cases,  except  the  gate  width,  the  standard  deviations  are 
taken  to  be  3%  to  10%  of  the  nominal  values  for  initial  designs.  The  gate  width  standard 
deviation  is  estimated  to  be  the  gate  length  standard  deviation  times  the  number  of  gate 
fingers.  The  variances  are  listed  in  Table  5.1. 
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Buried  Channel  Device 

Variable 

Variance 

Ls 

iB9 

flow 

IBQSBIfSSlil 

KEiSfl 

Nlow 

Nhigh 

VdD 

Uniform  Channel  Device 

Variable 

Variance 

Ls 

K&Q£MfZESH 

tchnl 

Nehnl 

BBMSM 

Vgg 

mKBSuajuwm 

■aw 

Ion  Imj 

planted  Device 

Variable 

Variance 

Ls 

■H 

4.0  X  10  ^  (m™)* 

E 

D 

Vgg 

fff'TtfVMira 

■339 

Table  5.1:  Disturbance  covariance  matrices  for  ion  implanted  device,  buried  channel  device, 
and  uniform  device 


Zload 

50+j0 


A  number  of  other  conditions  were  specified  for  the  experiments.  Three  harmonics  were 
used  during  harmonic  balance  calculations.  The  simulated  circuit  was  reduced  as  shown  in 
Figure  5.1.  The  impedances  presented  to  the  gate  and  drain  of  the  MESFET,  as  previously 
indicated,  were  50  +  jO  fi  at  all  harmonics.  During  each  yield  optimization,  the  number 
of  sample  devices  for  each  yield  estimate  was  100.  The  number  of  sample  devices  in  the 
presented  yield  histograms  is  500. 

In  Figures  5.2,  5.3,  and  5.4  we  show  the  initial  and  optimized  yield  histograms  for  ion 
implanted,  buried  channel,  and  uniform  channel  devices  when  small-signal  transducer  gain 
is  used  as  the  acceptance  criterion.  The  three  initial  designs  all  have  similar  initial  gain 
distributions  with  means  varying  from  5.1  to  5.9  while  the  standard  deviations  range  from 
0.7  to  1.1.  The  buried  channel  device  exhibits  the  best  improvement  with  a  distribution 
mean  increasing  to  16.0.  All  the  optimized  distributions  exhibit  greater  spread  than  their 
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Small-signal  transducer  gain 

Figure  5.2:  Simulated  histograms  of  ion  implanted  device  optimized  for  yield  relative  to 
small-signal  transducer  gain.  A  -  initial  design.  B  -  optimized  design. 

corresponding  initial  distributions.  The  standard  deviations  of  the  optimized  designs  range 
from  1.4  to  2.3.  The  increased  spread  results  from  the  optimized  designs  being  in  a  region 
were  the  gain  is  more  sensitive  to  perturbations  in  the  design  parameters.  For  all  three  devices 
the  optimizer  changes  the  gate  width  the  most.  These  changes  improve  the  matching  between 
the  transistor  and  the  50  (I  circuit.  This  large  change  in  gate  width  would  not  necessarily 
be  expected  for  a  tuned  circuit.  The  buried  and  uniform  channel  devices  also  exhibit  large 
changes  in  their  biases.  For  both  devices,  Vqg  and  Vdd  increswe.  This  bias  shift  is  to  a  region 
of  higher  transconductance.  The  increased  transconductance  improves  the  smaU-signal  gain. 

Figures  5.5,  5.6,  and  5.7  show  the  initial  and  optimized  yield  histograms  for  ion  implauited, 
buried  channel,  and  uniform  channel  devices  when  output  power  at  1  dB  gain  compression 
is  used  as  the  acceptance  criterion.  The  initial  designs  again  have  similar  performance.  The 
initial  design  distribution  means  range  from  253  to  316  mW  while  the  standard  deviations 
vary  from  55  to  71  mW.  The  uniform  channel  device  shows  the  greatest  increase  in  distribu¬ 
tion  mean  to  2283  mW.  However,  the  distribution  spread  is  more  than  twice  that  of  the  two 
other  optimized  designs.  In  all  three  cases,  the  optimized  designs  show  marked  improvement 
over  the  initial  designs.  Again,  the  optimizer  changes  the  gate  widths  of  all  three  devices  the 
most.  However  in  this  case,  the  increased  gate  widths  allow  more  current  to  flow  through 
the  devices  —  impedance  matching  becomes  a  secondary  consideration.  For  the  buried  and 
uniform  channel  devices  Vgg  and  Vdd  increase  considerably.  Vdd  is  also  markedly  increased 
for  the  ion  implanted  device.  The  biases  change,  along  with  the  channel  doping  and  gate 
width,  so  as  to  maximize  the  intersection  of  the  50  fl  load  line  and  the  device  I-V  curves. 

Yield  optimization  using  maximum  power-added  efficiency  as  the  acceptance  criterion 
was  also  performed.  Figures  5.8,  5.9,  and  5.10  show  the  initial  and  optimized  yield  his¬ 
tograms  for  ion  implanted,  buried  channel,  and  uniform  channel  devices.  The  initial  designs 
in  this  case  are  not  all  similar.  The  initial  uniform  channel  design  has  a  distribution  mean 
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Small-signal  transducer  gain 

Figure  5.3:  Simulated  histograms  of  buried  channel  device  optimized  for  yield  relative  to 
small-signal  transducer  gain.  A  -  initial  design.  B  -  optimized  design. 


Figure  5.4:  Simulated  histograms  of  uniform  channel  device  optimized  for  yield  relative  to 
small-signal  transducer  gain.  A  -  initial  design.  B  -  optimized  design. 
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Figure  5.5;  Simtilated  histograms  of  ion  implanted  device  optimized  for  yield  relative  to 
output  power  at  1  dB  gain  compression.  A  -  initial  design.  B  -  optimized  design. 
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Figure  5.6:  Simulated  histograms  of  buried  channel  device  optimized  for  yield  relative  to 
output  power  at  1  dB  gain  compression.  A  -  initial  design.  B  -  optimized  design. 
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Figure  5.7:  Simulated  histograms  of  uniform  channel  device  optimized  for  yield  relative  to 
output  power  at  1  dB  gain  compression.  A  -  initial  design.  B  -  optimized  design. 

of  19%,  whereas,  the  initial  ion  implanted  and  buried  channel  designs  have  distribution 
means  of  38.5%  and  40.0%,  respectively.  Yield  optimization  improved  the  uniform  channel 
design’s  distribution  mean  to  40.9%,  but  the  optimized  buried  channel  design  has  the  best 
performance  with  a  distribution  mean  of  47.6%.  We  note  no  improvement  in  the  ion  im¬ 
planted  design.  This  result  was  anticipated  since  the  ion  implanted  initial  design  is  based  on 
a  mature  industry  device  which  has  been  empirically  optimized  for  maximum  power-added 
efficiency.  As  was  the  case  with  the  gain  and  output  power  at  1  dB  gain  compression,  the 
largest  changes  occur  in  the  gate  width  and  the  DC  bias.  The  changes  in  gate  width  improve 
the  match  to  the  50  fi  circuit.  The  bias  point  shift  in  such  a  way  as  to  minimize  the  DC 
power  supplied  to  the  device  while  maximizing  the  intersection  of  the  50  fl  load  line  and  the 
device  I-V  curves. 

Different  optimum  device  designs  result  for  each  of  the  specified  performance  criteria. 
That  is,  an  optimum  PAE  design  is  different  from  an  optimum  design  for  either  maximum 
output  power  or  gain.  This  series  of  experiments  indicates  that  the  buried  channel  device  is 
the  best  device  structure  of  the  three  when  gain  and  power-added  efficiency  are  of  primary 
concern.  The  uniform  channel  device  gave  the  best  average  performance  for  output  power  at 
1  dB  gain  compression,  but  the  performance  distribution  exhibits  excessive  variance.  These 
result,  of  course,  only  apply  to  the  devices  embedded  in  a  50  D  circuit.  Different  results  are 
possible  when  circuit  tuning  conditions  are  considered. 
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Figure  5.8:  Simulated  histograms  of  ion  implanted  device  optimized  for  3deld  relative  to 
maximum  power-added  efficiency.  A  -  initial  design.  B  -  optimized  design. 
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Figure  5.9;  Simulated  histograms  of  buried  channel  device  optimized  for  yield  relative  to 
maximum  power-added  efficiency.  A  -  initial  design.  B  -  optimized  design. 
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Figure  5.10:  Simulated  histograms  of  uniform  channel  device  optimized  for  yield  relative  to 
maximum  power-added  efficiency.  A  -  initial  design.  B  -  optimized  design. 
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Chapter  6 
Conclusions 


The  NCSU  large-signal,  non-linear  MESFET  simulator  (TEFLON)  has  been  significantly 
advanced  and  enhanced  with  RF  sensitivity,  yield  determination,  and  yield  optimization 
routines.  The  physics  based  device  simulator  allows,  for  the  first  time,  determination  of  DC 
and  RF  performance  directly  in  terms  of  structural,  doping,  bias,  and  RF  circuit  tuning 
conditions.  Doping  profiles  can  be  determined  theoretically  using  the  Stanford  SUPREM 
3.5  process  simulator,  or  experimentally  using  a  physical  parameter  extraction  routine  that 
has  been  developed.  The  simulator  (TEFLON)  has  been  structured  to  link  directly  with 
SUPREM,  which  must  be  obtained  directly  from  Stanford.  The  physical  parameter  extrac¬ 
tion  routine  is  included  in  TEFLON. 

The  yield  simulator  uses  a  Monte  Carlo  method  that  makes  numerous  calls  to  the  device 
model,  until  an  accurate  estimate  of  the  yield  is  obtained.  Incorporation  of  an  advanced  near- 
global  optimization  algorithm  permits  the  initial  device  design  and  operating  conditions  to 
be  varied  until  a  maximized  yield  is  obtained.  In  this  manner  ’design  centering’  can  be 
efficiently  performed.  The  mathematical  optimizer  is  a  quasi-Newton  method  which  uses 
decreasing  gradient  scales  to  overcome  the  estimator’s  Monte  Carlo  noise.  The  optimizer 
permits  physical  parameters  to  be  adjusted  until  a  maximized  yield  in  any  specified  RF 
performance  measure  is  obtained.  Typically,  linear  gain,  RF  output  power  at  1  dB  gain 
compression,  and  maximum  power-added  efficiency  are  used  as  the  performance  measure. 

Use  of  the  simulator  for  yield  determination  and  design  centering  is  demonstrated.  It 
is  interesting  that  this  investigation  reveals  that  the  same  design  is  not  optimum  for  the 
different  performance  measures.  In  fact,  if  optimized  performance  with  maximum  yield  is  to 
be  obtained,  each  application  will  require  a  specialized  device  and  circuit  design. 

The  simulator  has  been  verified  to  be  quantitatively  accurate  for  simulation  of  MESFETs 
operating  at  C  and  X-bands.  Verification  tests  were  performed  using  device  data  supplied 
by  ITT,  Raytheon,  Texas  Instruments,  General  Electric  and  Hughes.  Excellent  agreement 
between  simulated  and  measured  data  was  obtained. 

An  early  version  of  the  simulator  and  User’s  Manual  was  offered  and  supplied  at  no 
cost  to  any  MIMIC  Program  participants  that  desired  a  copy.  Copies  were  delivered  to 
ITT,  Raytheon,  Texas  Instruments,  Hughes,  TRW,  Martin-Marietta,  Compact  Software, 
General  Electric,  Wright-Patterson  AFB,  and  Fort  Monmouth.  In  addition,  user  training 
was  provided  to  ITT  and  Raytheon,  also  at  no  cost.  Effective  use  of  the  simulator  requires 
some  learning,  but  the  effort  should  be  productive.  The  updated  version  of  the  complete 
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simulator  (TEFLON  5.0)  and  updated  User’s  Manual  are  available  from  NCSU, 

The  simulator  is  state-of-the-art  and  represents  a  significant  advancement  in  the  de¬ 
velopment  of  microwave  computer-aided  design  tools.  When  used  by  a  skilled  engineer  in 
industrial  applications,  the  simiilator  should  prove  useful  in  efforts  to  reduce  the  time  and 
costs  required  to  produce  advanced  MIMICs. 
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Appendix  A 

Yield  Estimate  Error 


The  true  yield  is 

1>= 

The  estimate  Y  has  the  correct  mean 

and  so  is  an  unbiased  estimate  of  <f).  The  variance  of  Y  is 

««r(r)  =  ((r  -  (y))’) 

is  equal  to 

var(y)  =  ^  -  ^))- 

But  the  disturbances  of  different  devices  are  uncorrelated  so  that 

’-<■-•(5')  =  E  E  -  ?)’)  =  ^ 

since  6ij  is  the  usual  Kronecker  delta  which  is  unity  when  i  =  j  and  is  zero  otherwise.  Now 

t;ar(0)={(0(*)-0)»)  =  (0*)-0* 
but  <l>{x)  is  either  unity  or  zero  so  that  =  <f>{xi)  and  so 

{<^(®)’}  =  (^(®))  =  1> 

and  therefore,  combining  these  results 
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Appendix  B 


Statistics  of  Disturbance  Vector 
Sequences 


We  now  show  that  our  Gaussian  generator  produces  disturbance  vector  sequences  with  the 
correct  statistics.  Let  y  =  x  +  v,  then  the  mean  of  y  is  x. 

{y)  =  (x),  +  Uw{z)^ 

because  the  expectation  operator  (.)  is  linear.  But  x  is  constant  so  that  (x)^  —  x.  By 
construction,  the  random  vector  z  is  zero  mean  so  that  (z),  =  0. 

Also,  the  covariance  of  y  is  C: 

cov(y)  =  ((y  -  (y).)  (y  - 


which  is 

cov(y)  =  ((Uwz)  (Uwz)^),  =  (Uwzz^w^U^),. 

But  U  and  w  are  constant  with  respect  to  the  expectation  operator  and  (z2^)z  =  I  the 
identity  matrix  since  z  is  a  vector  of  independently  distributed,  (Gaussian)  random  variables 
of  zero  mean  and  unit  variance,  so 

cov(y)  =  Uw(zz^),w^U^  =  Uww^U^. 

Finally  w  is  the  diagonal  matrix  of  the  square  roots  of  W,  so  that  =  W  and 

cov(y)  =  UWU'^, 

which  is  simply  the  originally  specified  covariance  C,  since  U  and  W  were  obtained  by  solving 
the  eigenproblem  problem  C  =  UWU^. 

Therefore  the  formula  y  =  x  +  Uwz  converts  a  vector  of  zero-mean,  unit-variance  inde¬ 
pendently  random  numbers  into  a  vector  of  process  disturbances  that  agree  with  the  specified 
mean  and  covariance  of  those  disturbances. 
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Appendix  C 

Sensitivity  Analysis  of  Ion  Implanted 
Device 


The  data  presented  in  this  Appendix  is  an  sensitivity  analysis  of  an  ion  implanted  MESFET 
device.  The  study  was  undertaken  to  identify  the  parameters  to  which  certain  performance 
measures  are  n;iost  sensitive.  The  performance  measures  considered  were  the  small-signal 
transducer  gain,  the  output  power  at  1  dB  transducer  gain  compression,  and  the  maximum 
transducer  power-added  efficiency.  These  measures  are  most  sensitive  to  gate  length,  gate 
width,  implant  dose,  implant  energy,  gate-source  bias  voltage,  and  drain-source  bias  voltage. 
During  the  simulations  for  which  the  DC  bias  voltages  were  not  analysis  parameters,  the 
gate-source  bias  voltage  was  set  to  0.51^,  and  the  drain-source  bias  voltage  was  set  to 
O.bVdtU- 
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Figure  C.2:  The  effect  of  gate  length  varia- 
Figure  C.I:  The  effect  of  gate  width  varia-  ^mall-signal  transducer  gain,  (b) 

tion  on:  (a)  smaU-signal  transducer  gam,  (b)  j  jg  compression,  and 

output  power  at  1  dB  gain  compression,  and  power-added  efficiency. 

(c)  maximum  power-added  efficiency. 
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Implant  Dose  (x  10**12  ions/cm**2) 


(a) 


(b) 


(c) 


Figure  C.3:  The  effect  of  implant  dose  varia¬ 
tion  on:  (a)  small-signal  transducer  gain,  (b) 
output  power  at  1  dB  gain  compression,  and 
(c)  maximum  power-added  efficiency. 
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Implant  Energy  (keV) 

(c) 

Figure  C.4:  The  effect  of  implant  energy  vari¬ 
ation  on:  (a)  small-signal  transducer  gain,  (b) 
output  power  at  1  dB  gain  compression,  and 
(c)  maximum  power-added  efficiency. 
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Figure  C.5:  The  effect  of  implant  cap  thick-  Figure  C.6:  The  effect  of  implant  anneal  tem- 
ness  variation  on:  (a)  small-signal  transducer  perature  variation  on:  (a)  small-signal  trans¬ 
gain,  (b)  output  power  at  1  dB  gain  com-  ducer  gain,  (b)  output  power  at  1  dB  gain 
pression,  and  (c)  maximum  power-added  effi-  compression,  and  (c)  maximum  power-added 
ciency.  efficiency. 
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(b) 


(c) 

Figure  C.7:  The  effect  of  implant  anneal 
time  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 


Gate  Source  Voltage  (V) 
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(b) 


(c) 

Figure  C.8:  The  effect  of  DC  gate-source  bias 
voltage  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.9:  The  effect  of  DC  drain-source  Figure  C.IO:  The  effect  of  parasitic  gate  re¬ 
bias  voltage  variation  on:  (a)  small-signal  sistance  variation  on:  (a)  small-signal  trans- 
transducer  gain,  (b)  output  power  at  1  dB  ducer  gain,  (b)  output  power  at  1  dB  gain 
gain  compression,  and  (c)  maximum  power-  compression,  and  (c)  maximum  power-added 
added  efficiency.  efficiency. 
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Drain  Resistance  (ohm) 
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Figure  C.ll:  The  effect  of  parasitic  drain  re¬ 
sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.12:  The  effect  of  parasitic  source  re¬ 
sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.13:  The  effect  of  parasitic  gate  in¬ 
ductance  variation  on:  (a)  smadl-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.14:  The  effect  of  parasitic  drain  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.15:  The  effect  of  parasitic  source  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 


Surface  Charge  Density  (cm**-2) 

(a) 


Surface  Charge  Density  (cm**-2) 

(b) 


Surface  Charge  Density  (cm**-2) 
(c) 


Figure  C.16:  The  effect  of  surface  charge  den¬ 
sity  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 
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Sensitivity  of  Maximum  PAE  Sensitivity  of  Pout  di  1  dB  GC  Sensitivity  of  Linear  Gain 
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Figure  C.17:  The  effect  of  density  of  implant 
cap  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power- added  efR- 
ciency. 
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Figure  C.18:  The  effect  of  electron  saturation 
velocity  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  C.19:  The  effect  of  electron  low-field 
mobility  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 


Figure  C.20:  The  effect  of  electron  critical 
electric  field  variation  on:  (a)  small-signd 
transducer  gadn,  (b)  output  power  at  1  dB 
gain  compression,  and  (c)  maximum  power- 
added  efficiency. 
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Figure  C.21:  The  effect  of  electron  high-field 
diffusion  coefficient  variation  on:  (a)  small- 
signal  transducer  gain,  (b)  output  power  at 
1  dB  gain  compression,  and  (c)  maximum 
power-added  efficiency. 
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Figure  C.22:  The  effect  of  gate  built-in  volt¬ 
age  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maLximum  power-added  effi¬ 
ciency. 
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Figure  C.23:  The  effect  of  Debye  length  varia¬ 
tion  on:  (a)  small-signal  transducer  gain,  (b) 
output  power  at  1  dB  gain  compression,  and 
(c)  maximum  power-added  efficiency. 
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Appendix  D 

Sensitivity  Analysis  of  Buried 
Channel  Device 


The  data  presented  in  this  Appendix  is  an  sensitivity  analysis  of  a  buried  channel  MESFET 
device.  The  study  was  undertaken  to  identify  the  parameters  to  which  certain  performance 
measures  are  most  sensitive.  The  performance  measures  considered  were  the  small-signal 
transducer  gain,  the  output  power  at  1  dB  transducer  gain  compression,  and  the  maximum 
transducer  power-added  efficiency.  These  measures  are  most  sensitive  to  gate  length,  gate 
width,  thicknesses  and  doping  densities  of  the  high  and  low  doped  regions,  gate-source 
voltage,  and  drain-source  voltage.  During  the  simulations  for  which  the  DC  bias  voltages 
were  not  analysis  parameters,  the  gate-source  bias  voltage  was  set  to  O.BVpo,  and  the  drain- 
source  bias  voltage  was  set  to  0.5Vd«M- 
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Sensitivity  of  Maximum  PAE  Sensitivity  of  Pout  or  1  dB  GC  Sensitivity  of  Linear  Gain 
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Figure  D.3:  The  effect  of  low  doping  region 
thickness  variation  on;  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  D.4:  The  effect  of  low  doping  den¬ 
sity  variation  on;  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 
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Sensitivity  of  Maximum  PAE  Sensitivity  of  Pout  fji  ’  1  dB  GC  Sensitivity  of  Linear  Gain 
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Figure  D.5:  The  effect  of  high  doping  region 
thickness  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 


Figure  D.6:  The  effect  of  high  doping  den¬ 
sity  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 
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Sensitivity  of  Maximum  PAE  Sensitivity  of  Pout  („  1  dB  GC  Sensitivity  of  Linear  Gain 
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Figure  D. 7:  The  efFec.  of  DC  gate-source  bias  Figure  D.8:  The  effect  of  DC  drain-source 
voltage  variation  on;  (a)  small-signal  trans-  bias  voltage  variation  on:  (a)  small-signal 
ducer  gain,  (b)  output  power  at  1  dB  gain  transducer  gun,  (b)  output  power  at  1  dB 
compression,  and  (c)  maximum  power-added  gain  compression,  and  (c)  maximum  power- 
efficiency.  added  efficiency. 
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Figure  D.9:  The  effect  of  parasitic  gate  re-  Figure  D.IO:  The  effect  of  parasitic  drain  re¬ 
sistance  variation  on:  (a)  small-signal  trans-  sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain  ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added  compression,  and  (c)  maximum  power-added 
efficiency.  efficiency. 


80 
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(c) 
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Figure  D.ll:  The  effect  of  parasitic  source  re-  Figure  D.12:  The  effect  of  parasitic  gate  in- 
sistance  variation  on;  (a)  small-signal  trans-  ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain  ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added  compression,  and  (c)  maximum  power-added 
efficiency.  efficiency. 
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Figure  D.13;  The  effect  of  parasitic  drain  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  D.14;  The  effect  of  parasitic  source  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  D.15:  The  effect  of  surface  charge  den-  Figure  D.16:  The  effect  of  electron  saturation 
sity  variation  on:  (a)  small-signal  transducer  velocity  variation  on:  (a)  small-signal  trans- 
gain,  (b)  output  power  at  1  dB  gain  com-  ducer  gain,  (b)  output  power  at  1  dB  gain 
pression,  and  (c)  maximum  power-added  effi-  compression,  and  (c)  maximum  power-added 
ciency.  efficiency. 
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Figure  D.17:  The  effect  of  electron  low-field 
mobility  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  D.18:  The  effect  of  electron  critical 
electric  field  variation  on:  (a)  small-signal 
transducer  gain,  (b)  output  power  at  1  dB 
gain  compression,  and  (c)  maximum  power- 
added  efficiency. 
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Figure  D.19:  The  effect  of  electron  liigh-field 
diffusion  coefficient  variation  on;  (a)  small- 
signal  transducer  gain,  (b)  output  power  at 
1  dB  gain  compression,  and  (c)  maximum 
power-added  efficiency. 
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Figure  D.20:  The  effect  of  gate  built-in  volt¬ 
age  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 
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Appendix  E 

Sensitivity  Analysis  of  Uniform 
Channel  Device 


The  data  presented  in  this  Appendix  is  an  sensitivity  analysis  of  a  uniform  channel  MESFET 
device.  The  study  was  undertaken  to  identify  the  parameters  to  which  certain  performance 
measures  are  most  sensitive.  The  performance  measures  considered  were  the  small-signal 
transducer  gain,  the  output  power  at  1  dB  transducer  gain  compression,  and  the  maximum 
transducer  power-added  efficiency.  These  measures  are  most  sensitive  to  gate  length,  gate 
width,  channel  doping  density,  channel  thickness,  gate-source  voltage,  and  drain-source  volt¬ 
age.  During  the  simulations  for  which  the  DC  bias  voltages  were  not  analysis  parameters, 
the  gate-source  bias  voltage  was  set  to  0.51^,  and  the  drain-source  bias  voltage  was  set  to 
O.SKuu- 
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Figure  E.3:  The  effect  of  channel  thick-  Figure  E.4:  The  effect  of  channel  doping  den- 
ness  variation  on:  (a)  small-signal  transducer  sity  variation  on:  (a)  small-signal  transducer 
gain,  (b)  output  power  at  1  dB  gain  com-  gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi-  pression,  and  (c)  maximum  power-added  effi¬ 
ciency.  ciency. 


89 


Gate  Source  Voltage  (V) 

(a) 


Gate  Source  Voltage  (V) 

(b) 


Gate  Source  Voltage  (V) 
(c) 


Figure  E.5:  The  effect  of  DC  gate-source  bias 
voltage  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.6:  The  effect  of  DC  drain-source 
bias  voltage  variation  on:  (a)  small-signal 
transducer  gain,  (b)  output  power  at  1  dB 
gain  compression,  and  (c)  maximum  power- 
added  efficiency. 
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Sensitivity  of  Maximum  PAE  Sensitivity  of  Pout  Or  1  dB  GC  Sensitivity  of  Linear  Gain 
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Figure  E.7:  The  effect  of  parasitic  gate  re¬ 
sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.8:  The  effect  of  parasitic  drain  re¬ 
sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  m^Lximum  power-added 
efficiency. 
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Figure  E.9:  The  effect  of  parasitic  source  re¬ 
sistance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.’.O:  The  effect  of  parasitic  gale  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.ll:  The  effect  of  parasitic  drain  in¬ 
ductance  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.12:  The  effect  of  parasitic  source  in¬ 
ductance  variation  on;  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figur.  E.13:  Th.  effect  of  surf«:e  charge  dec 
sity  v^atioB  oo:  (a)  emall-sigaal  Iraasdacer 
gain,  (b)  output  power  at  1  dB  gain  com¬ 
pression,  and  (c)  maximum  power-added  effi¬ 
ciency. 
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Fjgure  E.14:  The  effect  of  electron  saturation 
velocity  variation  on:  (a)  small-signal  trans¬ 
ducer  gam,  (b)  output  power  at  1  dB  eain 
efficSn^cy*^"’  niaximum  power-added 
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Figure  E.15:  The  effect  of  electron  low-field 
mobility  variation  on:  (a)  small-signal  trans¬ 
ducer  gain,  (b)  output  power  at  1  dB  gain 
compression,  and  (c)  maximum  power-added 
efficiency. 
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Figure  E.16:  The  effect  of  electron  critical 
electric  field  variation  on:  (a)  small-signal 
transducer  gain,  (b)  output  power  at  1  dB 
gain  compression,  and  (c)  maximum  power- 
added  efficiency. 
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Figure  E.17:  The  effect  of  electron  high-field  Figure  E.18:  The  effect  of  gate  built-in  volt- 
diffusion  coefficient  variation  on:  (a)  small-  age  variation  on:  (a)  small-signal  transducer 
signal  transducer  gain,  (b)  output  power  at  gain,  (b)  output  power  at  1  dB  gain  com- 
1  dB  gain  compression,  and  (c)  maximum  pression,  and  (c)  maximum  power-added  effi- 
power-added  efficiency.  ciency. 
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Figure  E.19:  The  effect  of  Debye  length  varia¬ 
tion  on:  (a)  small-signal  transducer  gain,  (b) 
output  power  at  1  dB  gain  compression,  and 
(c)  maximum  power-added  efficiency. 


97 


